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1.0  INTRODUCTION  AND  ABSTRACT 

This  Is  the  final  report  on  Contract  No.  N00014-76-C-0452:  "Environ- 
mental Effects  on  Fatigue  Crack  Initiation."  The  purpose  of  the  research  was 
to  develop  models  which  describe  the  effect  of  humidity  on  microscopic  fatigue 
failure  processes  In  aluminum,  and  lead  to  a description  of  the  scatter  In 
fatigue  lifetime  which  results  from  the  statistical  nature  of  the  fatigue 
failure  process.  In  order  to  develop  a comprehensive  model,  spanning  all  the 
sequenclal  failure  mechanisms  of  Importance  for  at  least  one  alloy,  we 
restricted  the  scope  of  our  research  to  alloy  and  loading  parameters  as 
fol lows: 


1)  To  materials  for  which  principal  crack  nucleation  occurs  at 
or  near  the  the  surface  at  Intermetal  11c  particles  - such 
as  is  the  case  for  2000  and  7000  series  aluminum  alloys. 

2)  To  peak  hardness  or  overaged  alloys. 

3)  For  constant  cyclic  amplitude,  R * omin/0max  = -1»  and  w1th 
maximum  amplitude  less  than  the  alloy's  0.2%  yield 
strength. 

4)  To  sinusoidal  loading  with  frequencies  in  the  range 
1 - 5 Hz. 

5)  To  materials  relieved  of  macroscopic  (or  long  range) 
residual  surface  stress. 

Our  research  addressed  mechanisms  in  the  Initiation  stage  of  fatigue 
failure.  Here  we  use  initiation  to  refer  to  all  processes  prior  to  the 
formation  of  a terminal  crack  (of  sufficient  size)  whose  subsequent 
propagation  can  be  predicted  by  conventional  fracture  mechanics.  For  the 
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aluminum  alloys  studied  Initiation  Is  comprised  of  three  processes  not 
heretofore  described  by  fracture  mechanics: 

1)  Nucleatlon  - The  formation  of  a microcrack,  typically  at  an 
Intermetalllc  particle  site. 

2)  Early  Microcrack  Growth  - The  growth  of  cracks  with  lengths 
on  the  order  of  the  grain  size,  with  propagation  rate 
substantially  influences  by  the  presence  of  grain 
boundaries. 

3)  Coalescence  - The  linking  of  microcracks  to  form  a terminal 
crack. 

Analytical  models  of  processes  1 and  2 have  been  developed  and  are 
the  basis  for  a Monte-Carlo  simulation  of  fatigue  failure  devised  to  predict 
the  lifetime  time  properties  of  the  aluminum  alloys  of  interest.  In  the  com- 
puter simulation  we  nucleate  surface  microcracks  which  grow  according  to  rules 
derived  from  the  models.  The  variation  in  lifetime  from  specimen  to  specimen 
is  related  to  statistical  fluctuations  in  alloy  microstructure  from  specimen 
to  specimen.  The  statistical  effect  of  crack  coalescence  on  lifetime  enters 
the  simulation  from  variation  in  relative  locations  of  the  microcracks.  In 
Section  2.0  we  describe,  in  outline  form,  the  general  sequence  of  events 
observed  in  fatigue  failure  of  2000  and  7000  series  aluminum  alloys  of 
Interest,  and  the  role  of  humidity  in  the  failure  process.  In  Section  3.0, 
detailed  experimental  results  are  discussed  along  with  their  mechanistic 
implications.  Finally,  in  Section  4.0,  results  of  modeling  predictions  are 
compared  to  experimental  data. 
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2.0  OUTLINE  OF  FATIGUE  FAILURE  SEQUENCE  TO  BE  MODELED 

Nucleation 

Crack  nucleation,  for  cyclic  loading,  occurs  as  the  result  of 
accumulated  damage  due  to  microplastic  deformation  of  the  alloy.  Owing  to 
lack  of  bulk  constraint  of  deformation,  the  surface  Is  the  most  vulnerable 
place  for  crack  nucleation.  Exceptions,  wherein  subsurface  nucleation  Is 
possible  Include:  1)  presence  of  a compressive  residual  stress  at  the  surface; 
or  2)  presence  of  large  subsurface  discontinuities  (voids,  Inclusions)  which 
can  act  as  sites  for  crack  nucleation. 

The  rate  of  crack  nucleation  at  the  surface  Is  influenced  by  the  type 
and  size  of  discontinuities  at  or  near  the  surface  which  can  act  as  stress 
concentrators,  and  hence,  as  sites  for  nucleation.  In  the  alloys  of  interest, 
intermetallic  particles  in  size  range  5 - 40  um  are  common  nucleation  sites. 
Occasional  nucleation  of  cracks  at  grain  boundaries  is  also  observed,  however, 
and  it  is  known  that,  even  in  the  absence  of  any  initial  surface  discontinui- 
ties, cracking  can  ultimately  begin  as  the  result  of  an  intrusion/extrusion 
type  mechanism.  It  is  a weakest  link  criteria  which  governs  which  of  the 
several  modes  of  crack  nucleation  will  dominate  [Rl]. 

For  a given  surface  discontinuity,  the  number  of  fatigue  cycles 
required  for  nucleation  is  determined  by  achievement  of  a critical  strain 
energy  density  in  the  region  of  the  discontinuity.  Cycles  to  nucleation  are 
reduced  if  the  discontinuity  size  is  large,  if  the  discontinuity  is  found 
within  a large  grain  permitting  large  slip  distances,  and  if  the  crystallo- 
graphic orientation  of  the  grain  is  favorable  [R2].  During  fatigue,  ductility 
and  hence  propensity  to  further  microplastic  deformation  of  the  alloy  surface 
layer  changes.  For  the  peak  hardness  alloys  studied,  the  trend,  due  to 
cutting  of  precipitates,  is  for  the  surface  to  soften.  In  the  presence  of 
humidity,  however,  surface  ductility  is  decreased  over  that  found  for  dry  air. 
As  a consequence,  one  observes  a reduced  rate  of  crack  nucleation  at 
intermetallic  sites  at  elevated  humidity.  Humidity  can  also  directly  alter 
the  properties  of  the  discontinuity  itself  leading  to  further  complexity.  In 

3 


I 


Rockwell  International 

Science  Center 


SC5050.1FR 

A1  2219-T851,  for  instance,  an  observed  effect  ">f  humidity  is  to  embrittle  the 
interface  between  intermetal  1 ic  and  matrix  [Rl].  An  additional  complexity 
ignored  in  the  present  research  is  the  effect  of  microscopic  residual  strain 
fields  associated  with  disconti nuites,  which  are  present  in  the  as  received 
material . 

In  Section  3.0  we  discuss  the  final  status  of  our  modeling  of 
nucleation  phenomena.  This  includes:  (1)  test  of  an  analytical  relationship 
between  minimum  cycles  to  nucleation  and  the  intermetallic  size,  the  size  of 
the  grain  containing  the  intermetallic,  and  the  cyclic  stress  amplitude;  (2)  a 
method  for  measuring  the  effect  of  humidity  on  ductility  of  the  surface  layer 
of  an  alloy,  and  (3)  observations  on  interrelationships  among  intermetallic 
composition,  structure,  relative  humidity  and  the  propensity  to  nuc’eate 
cracks  at  intermetallic  particle  sites. 

Early  Microcrack  Growth 

The  rate  of  propagation  of  surface  microcracks  is  substantially 
smaller,  on  an  average,  than  expected  from  conventional  fracture  mechanics. 

We  find  that  the  retarded  growth  is  principally  the  result  of  large  positive 
crack  closure  stresses  developed  by  the  microcracks.  Here  we  use  closure 
stress  in  a generic  sense  to  refer  to  all  processes  which  lead  to  a reduced 
cyclic  plastic  deformation  per  cycle  at  the  crack  tip.  We  relate  the  rate  of 
growth  of  the  microcracks  to  an  effective  stress  intensity  range,  AKeff,  where 
oc  (omax  - occ),  (i.e.,  the  cyclic  stress  range,  the  maximum  stress  minus 
the  closure  stress). 

The  closure  stress,  occ,  for  a given  crack  is  a function  of  location 
along  the  crack  front  and  changes  as  the  crack  grows,  leading  to  fluctuation: 
in  growth  rate.  Three  independent  phenomena  lead  to  positive  values  of  acc, 
and  it  is  the  mechanism  which  produces  the  largest  acc  at  a given  point  on  the 
crack  front  which  determines  occ  at  that  point  and  that  instant.  The  three 
mechanisms  are: 
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1)  Contact  of  the  opposing  fracture  surface  In  the  region  of 
the  intermetal  lie  particle.  The  contact  is  apparently  the 
result  of  a change  In  shape  of  the  fractured  Intermetal  11c 
particle.  occ  resulting  from  this  mechanism  Is  independent 
of  the  maximum  cyclic  amplitude,  omax.  Closure  from  the 
mechanism  Is  of  greatest  Importance  when  the  microcrack  is 
of  a size  close  to  the  Intermetal  lie  from  which  has 
nucleated.  As  the  crack  grows  occ  <*  (2c^_1,  where  2c  is 
the  crack  length. 

2)  Contact  of  opposing  fracture  surfaces  due  to  fracture 
surface  roughness.  We  have  not  studied  this  mechanism  in 
detail,  but  find  that  it  sets  a lower  limit  on  ucc  of  =&0.2 

°max* 

3)  Residual  tensile  strain  at  the  crack  tip.  The  tensile 
strain  develops  as  the  result  of  an  incomplete  reversal  of 
plastic  deformation  at  the  crack  tip.  The  mechanism 
dominates  the  closure  stress  phenomena  for  all  but  the 
earliest  stage  of  growth  of  microcracks. 

An  important  difference  between  microcrack  growth  and  crack  growth  in 
fracture  mechanics  specimens  is  that  for  microcracks  the  plastic  zone  size  is 
not  as  predicted  by  continuum  analysis.  Instead,  the  plastic  zone  size  is 
principally  limited  by  the  distance  between  a microcrack  tip  and  the  next 
neighboring  grain  boundary.  By  a series  of  experiments,  we  show  in 
Section  3.0  that  the  grain  size  at  the  crack  tip  can  be  used  to  predict  occ, 
arising  from  residual  strain  at  the  crack  tip,  and  hence,  microcrack 
propagation  rate.  We  show  that  a direct  effect  of  humidity  is  to  reduce  occ 
and  so  lead  to  increased  crack  propagation  rate. 

An  additional  complication  should  also  be  noted.  While  crack  closure 
stresses  are  the  prime  factor  in  microcrack  growth,  closure  stress  affects 
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only  Mode  I (i.e.,  tensile  opening)  crack  propagation.  As  a consequence.  If 
acc  ^ °max  ^ *s  not  uncommon  to  observe  Mode  II  propagation.  This  is  more 
common  in  the  grain  of  Initiation  where  the  crack  is  under  the  influence  of  a 
large  closure  stress  from  the  intermetall 1c  initiation  site.  It  Is  also 
observed  as  a microcrack  begins  to  propagate  Into  an  exceptionally  large 
grain,  for  which  the  residual  tensile  strain  developed  at  the  crack  tip  is 
large.  Propagation  typically  reverts  to  Mode  I as  the  closure  stress  drops. 
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3.0  SUMMARY  OF  EXPERIMENTS  AND  MODELS  CONCERNING  CRACK  NUCLEATION 
AND  EARLY  MICROCRACK  PROPAGATION 


Nucleation 

The  common  site  for  nucleation  of  cracks  at  an  Intermetallic  particle 
is  at  the  interface  between  the  matrix  and  the  Intermetallic  proper.  Nuclea- 
tion Is  observed  at  the  circumference  of  spherical  intermetal  lies  (Fig.  I)  and 
also  Interior  to  particles  containing  lamina  of  matrix  material,  such  as  In 
the  example  in  Fig.  2 for  A1  2219-T851.  In  the  latter  case,  cracking  commonly 
begins  within  the  particle  at  the  matrix-lamina  Interface.  In  the  2000  and 
7000  series  aluminum  alloys  the  laminated  internal  structure  of  the  inter- 
metallic particles  Is  found  In  association  with  particles  which  have  been 
present  during  rolling  (i.e.  which  are  elongated  In  the  rollim  direction). 

The  laminae  lie  approximately  perpendicular  to  the  rolling  direction  (R.D. ) 
and  fracture  at  the  lamina- Intermetal  11c  Interface  occurs  only  for  loading  in 
the  R.D.  As  a consequence,  low  Fe  content  A1  2219-T851,  which  contains  few 
spherical  Intermetal  lies,  has  a substantially  longer  fatigue  life  when  the 
loading  axis  is  In  the  long  transverse  direction.  Crack  nucleation  leading  to 
failure  occurs  at  grain  boundaries.  For  A1  2219-T851  the  intermetallic 
particle  structure  is  related  to  particle  composition,  which  also  affects  the 
sensitivity  of  the  particles  to  crack  nucleation  as  a function  of  humidity. 

Our  observations  on  the  Interrelationships  among  particle  composition, 
structure,  humidity,  and  propensity  to  nucleation  of  microcracks  have  been 
reported  In  detail  elsewhere.[Dl] 

A theoretical  expression  to  relate  the  minimum  number  of  fatigue 
cycles  to  crack  nucleation  at  Intermetallic  particles  to  the  alloy 
microstructure  In  the  vicinity  of  a particle  has  been  developed  with  support 
of  Rockwell  IRSD  funds. [R2]  The  main  result  of  the  analysis  Is  that  the 
minimum  number  of  cycles  to  nucleation,  N,  Is  given  by 
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Fig.  1 Circumferential  crack  at  0 phase  particle  in  A1  22 19-T85 1 . 
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C0  and  o0  are  material  constants.  As  Illustrated  by  Fig.  3,  W is  the  width  of 
the  intermetallic  normal  to  the  principal  stress  axis,  D is  the  maximum 
distance  of  slip  along  the  surface  at  approximately  45“  to  the  principal 
stress.  The  model  distinguishes  two  steps  to  crack  nucleatlon;  (a)  formation 
of  a crack  Inside  to  the  intermetallic;  (b)  propagation  into  the  matrix.  In 
Instances  of  surface  pitting  by  corrosion,  step  (a)  presumably  Is  not  re- 
quired, and  the  analysis  suggests  that  In  any  case  step  (b)  is  rate  control- 
ling for  all  but  the  lowest  cyclic  stress.  Equation  1 describes  step  (b).  As 
the  effect  of  crystallographic  orientation  has  been  Ignored  In  the  analysis, 
Eqn.  (1)  Is  expected  to  be  a lower  bound  of  the  number  of  cycles  to  nuclea- 
tlon, with  larger  N occurring  for  less  favorable  orientations.  In  the  future 
a more  complete  theory  must  also  consider  the  effect  of  triaxial  constraints 
on  N. 

We  have  tested  Eqn.  (1)  with  ONR  support  by  comparison  to  experi- 
mental results  obtained  for  Al  2219-T85I.[D2]  Plotted  in  Fig.  4 are  W vs  D 
values  for  nucleatlon  at  Intermetallic  particles  approximately  occurring  at 
three  Increments  in  fatigue  cycles,  N,  as  indicated.  Polished  flexural 
fatigue  specimens  of  Al  2219  utilized  were  given  a light  chemical  etch  so  that 
the  grain  boundaries  could  be  observed  in  a scanning  electron  microscope 
( SEM) . Dashed  lines  in  Fig.  4a-c  are  a fit  of  Eqn.  (1)  to  a single  W,  D 

value,  (x),  in  the  upper  right  corner  of  4a.  As  expected,  Eqn.  (1)  defines  a 

lower  bound  of  the  data.  The  scatter  towards  larger  W,  D values  is  attributed 

principally  to  variation  In  crystallographic  orientation  from  grain  to  grain. 

While  these  are  preliminary  results,  they  do  suggest  the  utility  of 
Eqn.  (1)  in  describing  the  effect  of  discontinuity  size,  and  grain  size  on 
cycles  to  crack  nucleatlon.  We  also  know  that,  at  least  qualitatively,  the 
stress  dependence  of  the  expression  is  good.  Before  similar  experiments  are 
done  to  rigorously  test  the  expression,  however,  the  desire  is  to  incorporate 
the  effect  of  crystallographic  orientation  into  the  expression.  Currently,  in 
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Fig.  3 Inter-metallic  width  and  slip  distance  parameters  important 
to  crack  nucleation. 
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Fig.  4 Range  of  intermetall ic  widths  and  slip  distances  leading  to 
nucleation  for  three  increments  in  fatigue  for  A1  2219-T851. 
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research  supported  by  Rockwell,  we  are  approaching  this  with  a Monte  Carlo 
simulation,  by  replacing  omax  with  an  effective  surface  stress  which  is  a 
function  of  the  orientation  of  the  <111>  type  slip  planes  in  a grain  relative 
to  the  external  stress  axis.  Computer  generated  scatter  plots  similar  to 
those  of  Fig.  4 have  already  been  obtained  with  the  technique. 

Early  Microcrack  Growth 

The  first  Indirect  evidence  of  the  large  closure  stresses  developed 
by  microcracks  was  the  observation  that  the  surface  microcracks  in  A1  2048  and 
A1  2219  did  not  completely  close  upon  relaxation  of  external  load.  [03,4],  It 
was  also  noted  that  the  residual  crack  openings  were  particularly  large  if  the 
cracks  were  of  grain  size,  suggesting  that  the  closure  stress  increased  as 
crack  lengths  reach  the  grain  size  - consistent  with  observed  retardation  in 
cracking  rate.  A method  was  developed  to  determine  the  “average"  closure 
stress  of  a microcrack  by  a compliance  measurement. [D4]  Using  a miniature 
flexural  fatigue  specimen  (Fig.  5)  the  opening  of  microcracks  was  measured  as 
a function  of  surface  stress  giving  results  as  those  illustrated  in  Fig.  6. 
"Closure"  defined  by  Fig.  6 is  considered  to  be  an  average  for  the  crack  in 
that  it  is  known  that  the  true  closure  stress  is  a function  of  location  along 
the  crack  front.  Hereafter  we  refer  to  the  average  closure  stress  determined 
by  such  a measurement  as  cfcc. 

In  a search  for  the  fundamental  mechanisms  of  closure  for  microcracks 
a series  of  experiments  were  done  to  relate  a'cc  to  the  shape  of  microcracks 
and  to  crack  opening  parameters  as  a function  of  crack  length.  Two  regions  of 
crack  closure  behavior  were  distinguished: 

a)  Short  cracks  for  which  crack  length  was  less  than  a factor 
of  two  larger  than  the  width  of  the  intermetal  lie  from 
which  nucleation  took  place.  In  this  circumstance,  Fcc  was 
found  to  be  independent  of  omax.  Crack  closure  resulted 
from  a misfit  of  fracture  surfaces  in  the  region  of  the 
intermetal lie. 
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Schematic  illustration  of  jig  used  for  loading  of  a 
miniature  flexure  specimen  in  an  SEM. 
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Longer  microcracks  for  which 


°cc  = C' 


C'  Is  a material  constant  and  6(o)/6(amax)  is  the  ratio  of  opening  at  the 
crack  center  at  zero  and  maximum  stress,  respectively.  Closure  in  this  case 
is  attributed  to  residual  tensile  strain  at  the  crack  tip.  For  a given  crack, 
6 ( o ) / 6 ( x ) is  a constant  and  hence  occ  is  a fixed  fraction  of  omax.  The 
measurements  leading  to  these  conclusions  have  been  reported  In  Met. 

Trans. [D4-7] 

While  determination  of  'occ  sufficed  to  define  mechanisms  and  trends 
in  closure  development  with  microcrack  growth,  It  was  necessary  to  develop 
quantitative  procedures  to  measure  the  true  closure  stress  as  a function  of 
location  along  the  crack  front,  as  a means  to  define  and  test  quantitative 
models  of  microcrack  growth.  Three  series  of  experiments  are  described  in 
3. 1-3.3. 

a)  It  was  shown  that  the  crack  closure  stress,  occ,  at  a 
surface  crack  tip  arising  from  residual  tensile  strain 
could  be  determined  quantitatively  from  measurement  of 
crack  tip  opening  displacement  (CT0D).[D8] 


b)  Comparison  between  rate  of  microcrack  growth  dc/dN,  and  an 
effective  stress  intensity  range,  4Ke^f<c(omax  - occ)  was 


made  for  A1  2219-T851 . 


acc  was  determined  at  increments  in 


crack  growth  by  CTOD  measurements.  We  found  that 


AUKeff)r 


gave  a good  prediction  of  the  microcrack  growth  rate,  where 
A and  m are  material  constants. [09] 
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c)  An  essential  element  of  the  deformation  process  for 
microcracks  Is  that  the  plastic  zone  size  Is  not  as 
predicted  from  continuum  mechanics.  Instead  the  zone  of 
significant  deformation  Is  determined  by  grain  boundary 
constraints  and  for  simple  grain  shapes  can  be  predicted 
from  the  distance  of  the  crack  tip  to  the  next  grain 
boundary  (Fig.  7).  Thus 


°cc  = a £ °max  * (4) 

where  a is  a material  constant  which  Is  sensitive  to  hu- 
midity. [DIO]  First  evidence  for  the  non-continuum  nature 
of  the  crack  tip  deformation  was  published  in  Met. 

Trans. [D7]  In  3. 1-3.2  we  expand  upon  the  discussion  of  the 
phenomena,  concluding  in  3.3  a comparison  between  CTOD  and 
distance  of  crack  tip  to  grain  boundary  which  is  the  basis 
for  Eqn.  (4). 

3.1  RELATIONSHIP  OF  SURFACE  MICROCRACK  TIP  OPENING  DISPLACEMENT 
TO  CRACK  CLOSURE  STRESS  FOR  A1  2219-T851  [DS] 

The  crack  tip  opening  displacement  (CTOD)  of  surface  microcracks,  60 
to  150um  in  length,  developed  during  fully  reversed  fatigue  loading  of 
A1  2219-T851  was  measured  and  compared  to  an  average  value  of  closure  stress 
for  the  cracks.  Scanning  electron  microscopy  was  used  to  measure  CTOD  at  the 
surface  tips  of  the  cracks  by  using  a flexurally  loaded  specimen  in  the  micro- 
scope. Corresponding  values  of  closure  stress  were  determined  by  an  empirical 
procedure,  based  upon  measurement  of  two  values  of  crack  opening  displacement 
at  the  crack  center.  The  average  closure  stress  of  a microcrack  was  found  to 
be  largest  for  cracks  having  largest  values  of  CTOD.  For  microcracks  having  a 
small  average  value  of  closure  stress,  the  CTOD  was  approximately  that 
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Fig.  7 Crack  parameters  of  importance  to  microcrack  growth. 
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expected  from  an  elastic  analysis  of  crack  opening  displacement  by  Green  and 
Sneddon. 

3.1.1  INTRODUCTION 

During  fatigue  crack  propagation,  large  cyclic  plastic  strains  occur 
within  a segment  of  the  plastic  zone  just  adjacent  to  the  crack  tip.  The 
existence  of  this  zone  was  predicted  by  R1ce,[R3]  and  has  subsequently  been 
confirmed.[R4,5]  More  recently,  the  distribution  of  strain  within  the  cyclic 
deformation  zones  of  several  alloys  has  been  measured  by  backscatter  electron 
channeling  by  Davidson  and  Lankford. [R6, 7]  Rice  also  noted  that,  on  an  un- 
loading cycle,  the  tensile  strain  from  the  prior  loading  cycle  was  not 
completely  reversed  by  compression,  and  resulted  in  crack  closure  at  a 
positive  value  of  externally  applied  load.  As  a consequence,  the  cyclic 
strain  per  cycle  in  an  element  of  material  ahead  of  the  crack  tip  is 
reduced.  Computer  assisted  modeling  of  crack  closure  by  Newman  has  been  based 
on  this  picture,  and  has  successfully  predicted  the  effect  of  deformation 
history  on  fatigue  crack  propagation  rate.[R8] 

Environmental  factors  also  affect  cracking  rate;  these  have  been 
studied  by  measurement  of  the  crack  closure  stress,  occ.[R9-ll]  Environment 
can  influence  the  ductility  of  the  material  at  the  crack  tip,  alter  the 
residual  tensile  strain  at  the  crack  tip  and,  hence,  change  the  amount  of 
cyclic  strain  experienced  deeper  in  the  cyclic  deformation  zone  where  the 
environment  may  not  penetrate.  If  the  most  important  effect  of  the 
environment  is  on  the  cyclic  strain  per  cycle,  it  is  reasonable  to  relate  the 
environmental  effect  on  cracking  rate  to  an  effective  intensity  range,  A«eff, 
at  the  crack  tip.[R12]  Akeff  a (omax-acc),  where  amax  is  the  maximum  cyclic 
stress.  In  essence,  Akeff  is  intended  to  describe  processes  that  alter  only 
the  cyclic  deformation,  although  this  is  recognizably  a simplification  since, 
for  instance,  crack  opening  and  closing  stresses  can  differ. [R13]  Modeling  of 
environmental  effects  on  crack  propagation  in  several  alloys,  using  AKeff,  has 
met  with  variable  amounts  of  success.  The  studies  are  not  straightforward, 
since  measurement  of  occ  is  still  controversial  [R14]  and  is  complicated 
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by  the  fact  that  the  closure  stress  can  vary  greatly  from  one  location  on  a 
crack  front  to  another. [R15] 

The  modeling  of  the  microscopic  aspects  of  fatigue  crack  propagation, 
and  of  crack  closure  will  continue  to  Improve  as  more  Is  learned  about  the 
Interrelationship  among  the  magnitudes  and  spatial  distribution  of  cyclic 
deformation  near  the  crack  tip,  the  effect  of  environment  on  the  deformation, 
and  the  relationship  of  the  cyclic  deformation  to  experimentally  determined 
values  of  crack  closure  stress.  In  the  present  work,  these  Interrelationships 
are  studied  by  using  surface  microcracks  of  the  grain  size  developed  during 
the  fatigue  of  Al  2219-T851.  Crack  tip  opening  displacements  (CTOD)  are 
measured  for  surface  cracks  with  lengths  on  the  order  of  the  grain  size,  and 
the  CTOD  values  are  compared  to  crack  closure  stresses.  This  Is  done  for  two 
values  of  relative  humidity,  for  fatigue  in  air.  Use  of  microcracks  for  the 
stu<1y  has  the  advantage  that  large  variation  in  closure  stresses  are  en- 
countered from  crack  to  crack,  which  provides  a wide  range  of  values  of 
closure  and  CTOD  for  comparison. 

3.1.2  EXPERIMENTAL  PROCEDURES 

An  Al  2219-T851  alloy,  with  a composition  of  6.3  25%  Cu,  0.3  wt%  Fe, 
0.2  wt%  Mn,  balance  Al,  was  used  for  the  Investigation.  The  alloy  has  prev- 
iously been  used  for  crack  Initiation  studies  and  Its  microstructure  Is  des- 
cribed elsewhere.[D4]  Blanks  from  which  specimens  were  fabricated  were  taken 
from  the  center  of  a 2.5  cm  thick  plate.  To  minimize  residual  surface  stress, 
machining  cuts  were  made  with  a decreasing  depth  of  cut  per  pass  ending  with  a 
final  cut  that  removed  20pm  of  material.  The  specimens  were  subsequently 
polished  using  a 0.05pm  alumina  powder.  A small  tapered  flexural  fatigue 
specimen  design  (Fig.  5)  was  utilized  which  permitted  a specimen  to  be  loaded 
to  a desired  tensile  surface  stress  within  a scanning  electron  microscope 
(SEM)  for  measurement  of  the  crack  tip  opening  displacements  and  closure 
stresses  of  microcracks  at  the  surface.  The  specimen  was  loaded  In  the  SEM  by 
displacing  the  small  end  of  the  specimen  a measured  amount  using  the  jig. 
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Test  specimens  were  fatigued  at  5 Hz  In  laboratory  air,  with  either  5 
or  45%  relative  humidity  (RH),  and  were  then  transferred  to  the  microscope  for 
study.  Loading  during  fatigue  was  fully  reversed  with  a maximum  cyclic 
stress,  of  90%  of  the  yield  strength  of  the  alloy  with  oyfe1(1  = 

350  MN/m2.  Specimens  were  positioned  with  the  surface  stress  axis  parallel  to 
the  alloy's  rolling  direction.  For  this  stress  orientation,  the  principal 
crack  Initiation  occurs  at  Intermetalllc  particles  located  near  the  surface; 
subsequent  crack  propagation  into  the  matrix  Is,  typically,  transgranular  and 
noncrystal lographlc.  A total  of  5xl03  fatigue  cycles  were  applied  to  each 
specimen  In  order  to  develop  numerous  microcracks  on  the  surface  with  surface 
lengths  in  the  range  of  10  to  150um.  The  longer  ones  of  these  were  Investi- 
gated in  the  present  study. 

A.  Determination  of  the  Average  Crack  Closure  Stress 

Conceptually,  the  closure  stress  at  the  tips  of  a microcrack  at  a 
surface  can  be  found  by  measuring  the  crack  opening  displacement  near  the  tip 
vs.  applied  surface  stress.  A technique  of  this  type  would  be  the  optical 
equivalent  of  an  Elber  gauge. [R12]  Practically,  however,  the  displacements 
Involved  are  too  small  for  accurate  measurement  of  the  closure  stresses  at  the 
tip  locations.  Instead,  in  prior  research  an  "average"  value  of  closure 
stress,  occ,  for  a crack  has  been  obtained  by  studying  the  dependence  of  the 
crack  opening  displacement  at  the  center  of  the  crack  (at  the  surface)  on 
applied  surface  stress. [D4]  The  closure  stress  thus  obtained  is  an  average  In 
the  sense  that  It  gives  a single  value  of  closure  stress  for  a crack,  even 
though  the  actual  closure  stress  Is  a function  of  location  along  the  crack 
front.  Average  closure  stress  is  defined  In  terms  of  a crack  opening  dis- 
placement vs.  applied  surface  stress  dependence,  as  Illustrated  in  Fig.  6. 

An  empirical  procedure  for  estimation  of  Fcc,  for  a surface  micro- 
crack, has  also  been  developed;[D4-6]  It  requires  measurement  of  only  two 
values  of  crack  opening  displacement  for  two  values  of  surface  stress.  The 
procedure  was  used  in  the  present  research  and  greatly  reduced  the  work 
required  to  determine  acc  for  a given  microcrack.  Before  the  estimation 
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procedure  Is  explained,  the  notation  used  to  describe  the  crack  opening 
displacement  In  later  sections  Is  described  (Fig.  8). 

The  magnitude  of  crack  opening  displacement,  of  course.  Is  a function 
of  both  location  within  a crack  and  of  the  applied  surface  stress,  o.  In  this 
paper,  displacement  measurements  are  made  only  at  the  alloy  surface,  as  a 
function  of  location  along  the  length  of  the  crack.  The  z axis  Is  normal  to 
the  stress  axis  and  Is,  approximately,  parallel  to  the  crack  path.  6(z,o)  Is 
the  crack  opening  displacement  at  the  surface  along  the  crack,  and  Is  measured 
parallel  to  the  applied  stress  axis,  z ■ 0 at  the  center  of  a crack.  The 
actual  closure  stress,  occ($).  Is  a function  of  location,  $,  along  the  sub- 
surface crack  front  denoted  by  the  dashed  line  In  Fig.  7.  The  surface  crack 
tip  at  z * -c  corresponds  to  $ » 0,  and  z * +c  to  $ * it.  These  are  the  two 
points  on  the  crack  front  for  which  the  crack  tip  opening  displacement  Is 
measured.  Equation  (5)  relates  occ  to  the  displacement  data.[D6]  For 
A1  2219-T851 


occ  * 3.81  omax  6(0,0)/6(0, 0,^ x ) . (5) 

The  parameters  6(0,0)  and  <5(0,omax)  are  the  crack  opening  displacements  at  the 
center  of  the  crack  for  a surface  stress  of  0 and  a tensile  stress  of  o,^, 
respectively. 

B.  Measurement  of  the  Crack  Tip  Opening  Displacement 

Measurements  of  the  crack  tip  opening  displacement  for  each  crack 
studied  were  made  at  surface  crack  tips  (corresponding  to  z ■ c and  z * +c  In 
Fig.  8)  using  micrographs  of  the  crack,  such  as  the  example  In  Fig.  9.  With 
the  specimen  loaded  to  omax,  the  crack  opening  displacement  was  measured  as  a 
function  of  position,  z,  along  the  length  of  the  crack.  The  crack  opening  at 
each  position  was  then  normalized  with  respect  to  the  maximum  opening 
displacement  that  occurred  at  or  near  the  center  of  the  crack.  An  example  of 
the  dependence  of  normalized  crack  opening  displacement  along  the  crack 
length,  obtained  with  this  procedure,  Is  shown  In  Fig.  10.  Crack  opening  dls- 
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placement  Is  not  symmetrical  to  z - 0,  attesting  to  the  Irregularity  of  the 
crack  path  and  an  asymmetry  In  crack  depth. 

CTOD  values  were  obtained  by  extrapolating  the  displacement  data  to 
the  two  crack  tips,  as  Illustrated  for  one  crack  tip  In  Fig.  11.  Because  of 
the  Irregular  structure  of  the  crack,  some  ambiguity  frequently  existed  on  how 
far  back  from  the  crack  tip  to  begin  the  extrapolation.  The  procedure 
selected  to  minimize  this  uncertainty  was  to  look  for  a point  within  0.5  to  3% 
of  the  total  crack  length  from  a tip  for  a sudden  change  In  d6(z)/dz.  An 
extrapolation  was  then  made  by  using  3 to  4%  of  crack  length  beyond  this 
point.  The  displacement  6(z)  was  measured  parallel  to  the  stress  axis, 
between  points  on  opposite  sides  of  the  fracture  that  could  be  Identified  as 
mating.  An  example  of  this  technique  Is  Illustrated  In  Fig.  12.  Opening 
displacement  measurement  must  be  done  with  some  care,  since  chips  and  other 
artifacts  In  the  opposing  fracture  surfaces  can  lead  to  an  erroneous  view  of 
how  the  surfaces  mate.  One  must  also  realize  that  a crack  may  not  open 
precisely  parallel  to  the  applied  stress  axis.  A suitable  measurement  also 
cannot  be  made  if  a crack  tip  Is  branched.  These  restrictions  limit  the 
number  of  cracks  for  which  reliable  CTOD  measurements  can  be  made. 

Ratios  of  CTOD  to  opening  displacement  at  z = 0 were,  approximately, 
Independent  of  the  surface  stress,  provided  that  the  surface  stress  was  close 
to  the  average  closure  stress  or  somewhat  larger.  the  subsequent  sections, 
therefore,  the  notations  6(c)/s( 0)  and  6 ( -c )/6 (0 ) are  used  to  denote  the  CTOD 
values  obtained,  which  are  normalized  to  the  opening  displacement  at  the  crack 
center  and  taken  for  sufficiently  large  surface  stress  so  that  the  ratios  are 
essentially  stress  Independent. 

3.1.3  RESULTS 

In  general,  the  normalized  CTOD  values  of  <5(-c)/6(0)  and  6(c)/6(0) 
obtained  for  the  two  surface  tips  of  the  same  crack  were  not  Identical.  As  an 
example.  Fig.  13  shows  two  surface  tips  for  the  same  crack  having 
substantially  different  values  of  CTOD.  It  was  at  first  thought  that  these 
differences  might  result  from  the  crack  being  especially  shallow  near  the  tip. 
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Fig.  10  Normalized  crack  opening  displacement  along  the  length 
of  the  microcrack. 
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Fig.  11  Extrapolation  of  normalized  crack  opening  displacement 
to  the  crack  tip  to  determine  a normalized  CTOD. 
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Fig.  12  The  stress  axis  is  parallel  to  arrows  which  indicate  a set 
of  matching  features  on  opposing  sides  of  the  fracture  sur- 
face used  to  measure  crack  opening  displacement. 
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Fig.  13  Two  surface  crack  tips  of  the  same  microcrack  illustrating 
large  variation  in  CTOD  observed  even  for  the  same  crack. 
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for  tips  with  small  CTOD;  this,  however,  does  not  appear  to  be  the  case.  The 
distribution  in  sub-surface  microcrack  shapes  was  studied  by  deliberately 
pulling  fatigued  specimens  to  failure  in  tension.  No  special  crack  depth 
variations  were  found  which  could  explain  the  observed  variation  In  CTOD. 

A series  of  measurements  were  then  done  in  which  the  values  6(c)/6(0), 
6(-c)/6(0),  and  occ  were  determined  for  a selection  of  surface  cracks  in  spe- 
cimens fatigued  In  air  for  either  5 or  45%  RH.  Relationships  were  then  sought 
between  the  normalized  CTOD  values  of  each  crack  and  its  corresponding  average 
closure  stress.  The  best  correlation  between  the  average  closure  stress,  occ, 
and  the  crack  tip  opening  displacement  was  found  in  comparing  the  average  of 
the  two  normalized  CTOD  measurements  of  a crack  to  occ.  This  comparison  is 
shown  In  Fig.  14,  where  the  averaged  normalized  CTOD  Is  defined  to  be: 

6(0)  " 76Tot  mui  * <6) 

Notice,  that  the  figure  includes  data  points  for  cracks  in  specimens  fatigued 
In  both  5 and  45%  RH.  A trend  is  observed  for  cracks  with  the  largest  value 
of  5(c)/6(0)  to  have  the  largest  occ,  and  no  statistically  significant  dif- 
ference between  the  relationship  for  5 and  45%  RH  is  apparent.  The  overt 
effect  of  humidity  was  that  at  45%  RH,  secondary  cracking  was  often  observed 
near  the  crack  tips,  such  as  the  example  shown  in  Fig.  15,  attesting  to  a 
decreased  ductility  of  propagation  for  the  higher  humidity.  CTOD  was  not 
measured  for  these  crack  tips. 

3.1.4  DISCUSSION 

For  the  A1  2219-T851  alloy  studied,  fatigued  with  the  stress  axis  in 
the  rolling  direction,  the  principle  crack  initiation  sites  are  at  intermetal - 
lie  particles  at  the  alloy  surface  (such  as  the  example  in  Fig.  9).  Crack 
closure  of  surface  microcracks  in  the  alloy  is  thought  to  arise  from  two 
mechan1sms.[D7]  Just  after  nucleation,  when  a microcrack  is  of  a length  com- 
parable to  the  intermetallic  from  which  it  propagates,  a misfit  between 
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Fig.  14  Relationship  of  average  normalized  CTOD,  5(c)/S(0),  to 
average  crack  closure  stress,  a .Jo 

LC  max 
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opposing  fracture  surfaces  produces  a crack  closure  for  positive  applied 
stress.  Such  misfit  Is  of  major  Importance  only  for  short  microcracks,  and  as 
a crack  grows,  the  deformation  behavior  along  the  crack  front  becomes  the 
factor  most  likely  to  determine  closure.  The  present  study  was  limited  to 
cracks  with  lengths  of  about  the  width  of  the  grains  normal  to  the  stress 
axis.  This  was  done  to  make  It  more  likely  that  crack  closure  resulted  from 
residual  tensile  strain  along  the  crack  front.  Values  of  5(c)/6(0),  found  for 
microcracks  In  A1  2219-T851,  range  from  0.09  to  0.33  as  plotted  in  Fig.  14. 

The  value  of  0.09  is  approximately  that  which  one  would  expect  for  the  elastic 
opening  of  a crack.  Some  of  the  problems  In  Interpretation  of  these 
observations  are  discussed  below. 


A.  Elastic  Limit  of  CTOD 

The  smallest  values  found  for  6(c)/6(0)  correspond  approximately  to 
that  expected  for  an  elastic  opening  of  an  elliptically  shaped  crack,  a case 
analyzed  by  Green  and  Sneddon. [R16]  This  result  is  obtained  if  one  approxi- 
mates the  surface  displacement  of  a semi -elliptical  surface  crack  by 


6 ( Z ) = 5(0) 


(7) 


where  Eq.  (7)  Is  that  obtained  by  Green  and  Sneddon  along  the  major  axis  of  an 
elliptically  shaped  crack  loaded  normal  to  the  plane  of  the  crack.  An  extrap- 
olation procedure  is  then  used  to  find  the  elastic  CTOD,  5e(c)/6(0),  from  the 
6 ( z ) vs.  z dependence  near  the  crack  tip,  as  shown  in  Fig.  16.  Naturally, 
there  Is  some  ambiguity  as  to  where  to  begin  the  extrapolation.  But,  if  the 
same  criteria  as  followed  with  real  microcracks  is  used,  a value  6e(c)/6(0) 
0.09  is  obtained.  This  is  in  agreement  with  the  smallest  values  of  6(c)/6(0) 
found  for  A1  2219- T851. 
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PERCENTAGE  DISTANCE  FROM  CRACK  TIP 


Fig.  16  Extrapolation  of  Eq.  (7)  to  determine  an  elastic  CTOD. 
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B.  Plastic  Deformation  Limit  on  CTOD 

Rice  [R3]  has  suggested  that  the  maximum  crack  tip  opening 
displacement  s(c,omax)  should  approximately  be 

5(c,amax)  = ^ Rp  * (8) 

for  cyclic  loading,  where  ey  Is  the  strain  at  the  yield  strength  and  Rp  Is  the 
diameter  of  the  plastic  zone.  For  the  range  in  CTOD  values  encountered  in  the 
present  work,  the  corresponding  calculated  values  of  Rp  extend  from  35  to 
130pm.  In  contrast,  the  plastic  zone  diameter,  calculated  using  the  poly- 
crystalline yield  strength  for  a 60pm  long  surface  crack  at  0.9  oyieid,  is 
24pm.  The  crack  tip  opening  displacements  and  associated  apparent  size  of 
zone  of  deformation  at  the  crack  tip  are,  therefore,  larger  than  expected  from 
the  continuum  theory.  As  Inferred  from  Eqn.  (8)  however,  the  range  In  plastic 
zone  sizes  is  approximately  the  same  as  the  range  In  the  width  of  the  grains 
In  the  alloy  In  the  direction  of  crack  propagation. 

Recent  refinements  of  Eq.  (8)  by  McKeeking  [R17]  to  incorporate  the 
effect  of  strain  hardening  on  CTOD  provides  corrections  much  too  small, 
although  in  the  right  direction,  to  explain  the  large  values  of  CTOD 
observed.  The  form  of  Eq.  (8)  is  not  convenient  for  evaluation  of  the  CTOD 
for  small  Rp,  since  the  CTOD  value  predicted  does  not  converge  to  the  elastic 
CTOD.  If  the  elastic  evaluation  is  taken  as  a lower  bound,  and  Rp  from 
Eq.  (8)  as  an  upper  bound,  the  observations  of  this  study  are  consistent  with 
a range  in  Rp  values  from  0 to  130  pm. 

C.  Further  Considerations 

The  relationship  between  CTOD  and  closure  stress,  found  for  micro- 
cracks of  the  grain  size,  is  substantially  different  from  results  previously 
reported  for  macrocracks  in  fracture  mechanics  specimens.  Langford  and 
Davidson  [R8]  for  instance,  have  observed,  a decreasing  CTOD  at  maximum  stress 
as  a crack  was  propagated  into  a zone  producing  increased  closure.  Some  years 
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ago,  Adams  [R19]  pointed  out  that  this  variation  should  be  expected  If 
residual  tensile  strain  In  the  wake  of  the  crack  was  responsible  for  the  crack 
closure  phenomenon.  In  contrast,  for  cracks  of  the  grain  size,  the  trend 
found  here  Is  for  closure  stress  to  Increase  with  increased  CTOD.  Previous 
studies  of  the  closure  phenomena  have  found  other  ways  In  which  the  behavior 
of  microcracks  differs  from  that  for  macrocracks.  Some  of  the  prior 
observations  for  microcracks  were  that: 

1.  6 ( z = 0)  <x  o , for  o > Fcc;  whereas  from  continuum  theory 

6 ( z * 0)  Is  expected  to  be  noticeably  nonll nearly  related  to 
stress. [D7] 

2.  7CC  a omax;  where  from  continuum  theory  occ  a o^gX  is 
expected. [D7] 

3.  occ  can  be  estimated  from  the  geometry  of  the  crack  path  as  seen 
from  the  surface,  and  cracks  differing  In  size,  but  not  shape, 
tend  to  have  the  same  closure  stress. [04] 

4.  Tensile  spike  overloading  Increases  Fcc  In  proportion  to  the 
magnitude  of  the  overload,  with  no  measureable  crack 
propagation. [07] 

To  this  list  can  be  added  the  latest  results: 

5.  occ  Increases  as  the  average  normalized  CTOD  of  the  surface 
crack  tips  increase. 

6.  For  cracks  with  small  closure  stress,  the  average  normalized 
CTOD  of  the  surface  tips  Is  approximately  what  would  be  observed 
if  there  were  no  plastic  deformation  at  the  crack  tips. 
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7.  The  large  values  of  CTOO  observed  for  some  of  the  microcracks 
are  substantially  larger  than  expected  from  continuum  theory, 
and  plastic  zone  sizes  estimated  for  these  are  approximately  the 
size  of  the  largest  grains  In  the  material. 

An  analytic  theory  to  describe  these  observations  remains  to  be  ad- 
vanced. Qualitatively,  however,  substantial  consistency  between  the  obser- 
vations Is  achieved  if  it  Is  assumed  that  the  size  of  the  plastic  zone  (or  at 
least  the  zone  of  significant  deformation)  at  a surface  tip  of  a microcrack  Is 
approximately  the  distance  of  the  crack  tip  to  the  nearest  grain  boundary,  and 
that  the  zone  size  Is  consequently  Independent  of  the  applied  stress.  One  can 
only  speculate  as  to  how  this  circumstance  might  arise,  but  the  factors  unique 
to  the  microcracks,  studied  here,  are:  the  plastic  zone  size  expected  from 
continuum  theory  Is  on  the  order  of  the  grain  size;  because  of  the  small  crack 
size,  the  plane  stress  condition  at  the  surface  crack  tips  dominates  the  de- 
formation behavior;  and  achieving  even  a comparatively  small  stress  Intensity 
entails  application  of  a large  surface  stress  to  a specimen. 

To  explain  the  closure  observations  for  microcracks  accumulated  In 
this  and  In  prior  work,  the  following  model  Is  suggested.  Except  for  special 
circumstances  where  fracture  surface  roughness  leads  to  closure,  the 
Inflection  In  the  crack  opening  displacement  vs  applied  surface  stress 
dependence,  interpreted  as  closure  for  microcracks,  is  not  the  result  of  a 
physical  closure  of  the  crack  In  the  wake  of  the  crack  front.  Instead,  the 
departure  from  elastic  behavior  on  unloading  Is  the  result  of  a cessation  of 
reverse  yielding  In  the  plastic  zone.  (As  a parenthetical  note.  It  should 
still  be  useful  to  describe  crack  propagation  in  terms  of  Meff  for  this 
case.  The  so-called  "crack  closure"  has  the  same  effect  of  reducing  the 
cyclic  deformation  per  cycle  experienced  by  material  at  the  crack  tip.)  The 
crack  opening  displacement-applied  stress  dependence  typically  observed  Is 
approximately  linearly  elastic  for  o > occ,  because  the  size  of  the  plastic 
zone  Is  Independent  of  the  applied  stress.  One  also  expects  from  this  that 
the  CTOD  <x  a for  5 > occ  leading  to  a closure  stress  that  Is  proportional 
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to  omax,  as  is  observed.  Humidity  reduces  the  ductility  of  the  alloy  in  the 
grain  sized  plastic  zone,  leading  to  lower  CTOD,  a smaller  7CC,  and,  hence,  to 
accelerated  crack  propagation  rates.  Rapid  variation  in  crack  propagation 
rate,  as  a microcrack  grows,  are  principally  the  result  of  changing  closure 
stress  along  the  crack  front  as  the  crack  position  changes  relative  to  the 
grain  boundaries.  This  picture  is,  of  course,  more  complicated  if  crack 
branching  occurs.  A dimensionally  based  argument  can  be  made  that 


amax 


o 

7c~ 


(9) 


where  zQ  is  the  distance  of  a surface  crack  tip  to  the  closest  grain  boundary 
in  the  direction  of  crack  propagation,  and  2c  is  the  crack  length. 

Naturally,  there  is  a substantial  amount  of  speculation  in  the  model 
and  further  investigation  Is  required  to  test  it.  Two  opportunities  in  this 
regard  are:  (1)  compare  dc/dN  to  AKe^^,  by  measuring  the  microcrack  closure 
stress  using  the  CTOD  relationship  of  Fig.  14,  and  (2)  look  for  a relationship 
between  CTOD  and  the  distance  of  microcrack  tips  to  neighboring  grain 
boundaries.  Development  of  an  analytic  model  of  the  deformation  properties  of 
cracks  with  plastic  zones  of  fixed  size  is  also  ultimately  essential  for  com- 
plete Interpretation  of  the  observations  of  these  studies. 

An  additional  coment  is  required.  It  might  be  expected  that  if 

6(c)/5(0)  = 0.09,  then  occ  = 0.  Values  of  occ  = 0.2  are  more  typical,  as  per 

Fig.  14.  It  is  suggested  that  crack  closure  due  to  misfit  of  the  opposing 
fracture  surfaces  (l.e.,  fracture  surface  roughness)  places  a lower  bound  on 
occ.  The  presumption  Is  that  the  mechanism  which  gives  the  largest  closure 
stress  at  a point  on  a crack  front  which  determines  the  closure  stress  at  that 

point.  A value  of  0.2  for  Fcc  is  in  the  range  commonly  found  for  closure 

stresses  in  2000  series  aluminum  alloys  for  large  cracks  in  fracture  mechanics 
specimens. 
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3.1.5  SUMMARY 

Crack  tip  opening  displacement  (CTOD)  measurements  were  made  for  sur- 
face microcracks  with  lengths  of  aproxlmately  the  grain  size  developed  during 
the  fatigue  of  A1  2219-T851.  Correlations  were  sought  between  the  two  CTOD 
values  obtained  for  each  crack  for  the  surface  crack  tips,  and  the  average 
closure  stress  of  the  crack.  Average  closure  stress  was  obtained  by  using 
empirical  relationship  based  on  two  measurements  of  the  crack  opening  dis- 
placement at  the  crack  center  — one  at  zero  and  the  other  at  maximum  applied 
surface  stress.  The  average  closure  stress,  occ.  Is  an  average  In  the  sense 
that  a single  closure  stress  value  Is  obtained  by  the  measurement  procedure, 
even  though  It  Is  expected  that  the  actual  closure  stress  varies  with  location 
along  the  crack  front.  Cracks  with  the  largest  values  of  CTOD  exhibited  the 
largest  values  of  occ.  This  trend  was  strongest  when  a comparison  was  made 
between  acc  and  the  average  of  the  two  CTOD  values  obtained  for  a crack,  with 
CTOD  normalized  to  the  opening  In  the  crack  center.  A model  was  proposed  to 
explain  these  and  prior  results  for  microcracks,  results  which  differ  substan- 
tially from  closure  observations  for  large  cracks  In  fracture  mechanics  speci- 
mens. The  key  element  of  the  model  Is  the  suggestion  that  for  microcracks, 
the  plastic  zone  size  (or  at  least  the  zone  of  significant  deformation)  at  the 
surface  Is  approximately  given  by  distance  between  crack  tip  and  the  grain 
boundary  toward  which  a crack  Is  propagating.  Furthermore,  it  Is  noted  that  if 
the  closure  stress  at  the  surface  tip  Is  proportional  to  the  plastic  zone 
size,  much  of  the  closure  observations  unique  to  microcracks  are  qualita- 
tively explained.  An  imnediate  application  of  the  model  Is  to  permit  an 
experimental  comparison  between  the  rate  of  propagation  of  surface  microcracks 
during  fatigue,  and  an  effective  stress  Intensity  range,  The  closure 

stress  values  needed  to  complete  the  calculations  of  can  be  obtained  by 

CTOD  •measurements. 
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3.2  CRACK  CLOSURE  STRESS  EFFECTS  ON  THE  RATE  OF  PROPAGATION  OF  SURFACE 

MICROCRACKS  DURING  THE  FATIGUE  OF  A1  2219-T851  [09] 

The  rate  of  propagation  of  surface  microcracks  of  approximately  25- 
100  um  In  length,  developed  during  fatigue  loading  of  A1  2219-T851,  has  been 
measured  and  compared  to  predicted  cracking  rates  using  an  effective  stress 
Intensity  range,  analysis.  Crack  closure  stress  and  crack  shape  factor 

parameters  needed  to  evaluate  AKeff  for  individual  crack  tips  were  obtained 
from  measurement  of  crack  tip  opening  displacement  and  crack  compliance, 
respectively.  Although  the  aKeff  description  of  microcrack  growth  Is 
recognizably  simple,  It  Indicates  that  crack  closure  contributes  significantly 
the  rate  of  early  propagation  of  microcracks  In  the  alloy. 

3.2.1  INTRODUCTION 

As  a result  of  fatigue,  unusually  large  postlve  crack  closure 
stresses  are  developed  by  surface  cracks  of  lengths  comparable  to  the  grain 
size  In  aluminum  alloys. [D3-6]  While  a trend  towards  reduced  rate  of  propaga- 
tion of  such  microcracks  with  Increased  closure  has  been  noted,  a quantitative 
comparison  between  measured  closure  stress  and  cracking  rate  remains  to  be 
made.  As  with  prior  exerlments  for  crack  closure  propagation  in  fracture 
mechanics  specimens,  It  Is  desirable  to  determine  an  effective  stress  In- 
tensity range,  aKeff,  to  describe  the  cracking  process.  AKe^  is  defined  to 
be  proportional  to  the  cyclic  stress  range  experienced  at  the  crack  tip 
(l.e.  to  (anwx"°cc^’  or  the  max1mum  cyclic  stress  minus  the  crack  closure 
stress) .[R12]  A problem  heretofore  In  evaluating  AKeff  for  microcracks  is 
that  occ  Is  a function  of  location  along  the  crack  front.  Crack  closure 
values  previously  obtained  for  microcracks  have  been  averages,  termed  o"cc 
hereafter,  obtained  from  compliance  measurements  made  at  the  center  of  surface 
microcracks. [D4]  Such  a measurement  determines  a single  value  of  "closure 
stress"  for  a crack  (this  Is  the  sense  in  which  "average"  is  meant)  and  hence 
suffices  only  for  the  study  of  trends  In  closure  stress  development  with  crack 
propagation. 
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Recently,  a correlation  has  been  found  for  A1  2219-T851,  between  the 
crack  tip  opening  displacement  (CTOO)  of  a microcrack  tip  and  occ,  with  occ 
increasing  with  Increasing  CT0D.[D8]  In  the  present  research  AKeff  values  at 
surface  microcrack  tips  in  A1  2219-T851  are  obtained  with  closure  stresses 
estimated  by  measuring  CTOO.  It  has  been  suggested  that  observed  variations 
In  occ  and  CTOO  as  a microcrack  propagates,  reflects  changing  constraints  on 
the  plastic  zone  size,  or  at  least  the  zone  of  significant  deformation,  with 
distance  between  the  microcrack  front  and  grain  boundaries. [D8]  The  crack 
shape  factor  (depth/  length)  needed  to  complete  the  determination  of  Meff  Is 
found  by  measurement  of  the  compliance  of  the  microcrack.  Using  these 
techniques,  the  rate  of  crack  propagation  of  microcracks  developed  during 
fatigue  In  A1  2219-T851  Is  compared  to  AKeff,  as  the  local  closure  stresses 
change  with  crack  propagation. 

3.2.2  EXPERIMENTAL  PROCEDURE 

The  A1  2219-T85I  alloy  used  for  the  investigation  had  a composition 
of  6.3  wt%  Cu,  0.3  wt%  Mn,  0.3  wt%  Fe,  balance  Al.  Miniature  flexural  fatigue 
specimens  were  machined  from  the  center  of  a 2.5  cm  thick  rolled  plate.  To 
reduce  residual  surface  stress,  the  cutting  depth  during  final  machining  was 
decreased  with  each  pass,  ending  with  20ym  of  material  removed  per  cut. 
Finally,  specimens  were  polished  with  a Q.05um  alumina  powder.  Results  re- 
ported here  were  obtained  using  a flexural  fatigue  specimen  geometry  and  a 
loading  jig  (Fig.  5)  designed  for  use  in  a scanning  electron  microscope 
(SEM).  The  specimen  geometry  utilized  provides  a uniform  tensile  stress  over 
the  specimen  surface,  whose  magnitude  is  determined  by  measurement  of  dis- 
placement of  the  specimen  end.  Loading  of  the  specimen  In  the  microscope  is 
required  for  measurement  of  the  microcrack  tip  opening  displacements. 

Each  test  specimen  was  fatigued  (R  = -1)  at  5 Hz  in  air  at  room  tem- 
perature, and  at  a relative  humidity  of  40%,  with  amax  = 0.9  Oy-jg^.  Surface 
stress  calibration  in  the  specimen  gauge  section  was  obtained  by  strain 
gauges.  The  absolute  stress  accuracy  is  estimated  to  be  ±5%,  the  relative 
accuracy  from  one  cracking  area  to  another  in  on  the  order  of  ±1%.  Material 
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orientation  was  with  the  rolling  direction  parallel  to  the  principle  stress 
axis.  For  this  orientation,  the  alloy  yield  stress  was  350  MN/m^  and  the 
crack  nucleatlon  was  at  Intermetalllc  particles  at  or  near  the  surface, 
followed  by  transgranular  crack  propagation. 

Microcrack  propagation  data  were  taken  as  follows:  A specimen  was 
given  2.5xl03  fatigue  cycles  In  air  and  was  then  transferred  to  the  SEM. 
Microcracks  were  located,  each  having  a surface  length  of  approximately 
25um.  These  were  microcracks  for  which  propagation  into  the  matrix  was 
already  well  established. 

High  resolution  micrographs  were  taken,  along  the  length  of  each 
crack  studied,  with  the  specimen  under  load:  a)  concentrating  particularly  on 
the  two  crack  tips  for  use  In  estimation  of  the  crack  tip  closure  stress  from 
CTOD;  and  b)  around  the  center  of  the  crack  for  use  In  estimating  both  clo- 
sure stress  and  the  crack  shape  factor.  With  these  measurements  completed, 
the  specimen  was  returned  to  the  air  environment  for  additional  fatigue.  The 
sequence  of  cracking  parameter  measurements  and  fatigue  was  repeated  in  the 
same  manner  in  Increments  of  500  cycles  until  5xl03  cycles  had  been  applied. 

At  this  point,  the  cracks  being  observed  were  approximately  lOOpm  long.  The 
resulting  data  are  used  to  compare  measured  crack  propagation  rate,  dc/dN,  and 
aKgff.  Before  discussion  of  these  results.  It  is  advantageous  to  briefly  des- 
cribe the  calibration  procedures  utilized  In  estimating  the  crack  tip  closure 
stresses  and  crack  shape  factor. 

A.  CTOD  Estimate  of  Crack  Tip  Closure  Stress 

An  empirical  relationship  has  previously  been  noted  between  the  crack 
tip  opening  displacement  determined  under  load  and  crack  closure  stress  for 
microcracks  developed  by  the  fatigue  of  A1  2219-T851 .[D8]  While  theoretical 
arguments  can  be  advanced  for  the  form  of  a general  relationship  between 
closure  stress  and  CTOD,  the  functional  relationship  between  these  quantities 
has  been  demonstrated  for  an  R * -1,  omax  ■ 0.9  ay1eid  test  condition.  In  the 
current  experiment,  specimens  were  fatigued  for  this  same  condition. 
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The  nomenclature  used  to  describe  the  crack  opening  displacement 
along  a microcrack  is  illustrated  by  Fig.  8.  5(z,o)  is  the  opening  at  point  z 

along  the  crack  length,  at  surface  stress  o,  and  is  measured  parallel  to  the 
applied  stress  axes.  Mating  features  on  opposite  sides  of  the  crack  opening 
are  identified  and  used  to  define  the  crack  opening.  6(c,a)  is  the  CTOD  and 
is  obtained  by  extrapolation  of  6(z,a)  to  z = c. 

Equation  (10)  is  used  to  estimate  the  closure  stress  at  a crack  tip 
from  the  CTOD. [Taken  from  Fig.  14.] 

°cc  /°max  = 6(c,a)/5(0,o).  (10) 

A requirement  on  the  validity  of  Eqn.  (10)  is  that  o exceed  acc.  Therefore, 
it  is  most  convenient  to  use  a = a,^.  The  quantity  6(0, a)  is  the  opening, 
for  the  same  a,  at  the  crack  center.  It  should  also  be  recognized  that 

Eqn.  (10)  is  an  empirical  description  of  microcrack  closure  stress  spanning 

two  different  mechanisms  of  closure  for  microcracks:  1)  closure  due  to 
fracture  surface  roughness  of  importance  for  occ  ~ 0.2  amax;  and  2)  closure 
due  to  residual  tensile  strain  at  the  crack  tip,  important  for  occ  > 0.2 

amax*  when  the  tensi1e  residual  at  a crack  tip  is  small  6(c,a)/6(o,a)  ~ 0.1, 

corresponding  to  the  elastic  limit  of  CTOD.  In  this  case,  surface  roughness 
leads  to  a closure  stress  of  ~ 0.2  amax  for  A1  2219-T851. 

B.  Crack  Shape  Factor  Calibration 

To  complete  the  calculation  of  AKeff,  it  was  necessary  to  determine 
approximately,  the  subsurface  shape  of  an  observed  microcrack  as  well  as  the 
manner  in  which  the  shape  changed  with  increment  in  propagation.  It  is  ex- 
pected from  an  analysis  by  Green  and  Sneddon  [R16]  that  the  quantity 
6(0,  amax)/2c  will  approximately  be  a constant  for  a given  shape  factor, 
a/2c.  An  experimental  calibration  was  performed,  as  follows,  to  relate  the 
two  quantities.  A specimen  of  A1  2219-T851  containing  surface  cracks  induced 
by  fatigue  was  pulled  in  tension  to  failure.  Numerous  microcracking 
thumbnails  were  exposed  from  which  the  distribution  of  shape  factors  was 
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measured.  One  such  crack  Is  shown  In  Fig.  17.  There  tends  to  be  some 
distortion  In  the  apparent  shape  of  the  crack  due  to  tearing  In  the  region  of 
the  crack  tips  near  the  surface.  Consequently,  the  crack  shape  was  measured 
using  markers  of  crack  front  location  within  the  fracture  surface,  such  as  at 
the  arrow.  The  distribution  In  observed  crack  shapes  as  a function  of  crack 
length  Is  shown  In  Fig.  18. 

A relationship  between  a/2c  and  5(0,  amax)/2c  was  obtained  by 
comparing  Fig.  18  to  a plot  of  the  scatter  In  6(0,  Omaxl/Zc  with  2c  obtained 
for  A1  2219-T851  in  prior  research. [D5]  The  mean  value  of  a/2c  observed  was 
approximately  0.362  and  led  to  a linear  correction  to  a/2c  for  variation  In 
6(0,  omax)/2c  given  by  Eqn.  (11). 


= 0.362  + 25.1 

While,  owing  to  the  comparatively  small  variation  in  crack  shapes  around  the 
mean.  It  might  have  sufficed  to  use  the  mean  In  calculating  AKeff  for  all 
microcracks,  use  of  Eqn.  (11)  was  found  to  decrease  the  apparent  data  scatter 
In  comparison  of  dc/dN  to  AKeff. 

3.2.3  RESULTS 

The  requirement  for  accurate  measurement  of  CT00  substantially  limits 
the  number  of  cracks  suitable  for  stuqy.  A key  element  is  that  one  must  be 
able  to  Identify  a number  of  matching  features  on  opposite  sides  of  the  crack, 
for  definition  of  6(z,o)  near  the  crack  tip.  While  numerous  cracks  were 
observed,  data  shown  here  are  limited  to  four  selected  cracks  for  which 
reliable  5(z,o)  values  were  obtained  for  a major  portion  of  the  crack  propa- 
gation. Crack  #1  (Fig.  19)  is  Included  because  It  had  the  largest  observed 
range  In  dc/dN.  In  the  later  stage  of  fatigue,  the  lower  crack  tip  of  crack 
#1  branched  In  passage  of  an  intermetalllc  particle  (Fig.  20).  This  reduces 
the  accuracy  of  the  calculated  AKeff  for  that  tip  for  the  4.5xl03  and  5xl03 
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Fig.  17  Fatigue  crack  thumbnail  used  to  measure  crack  shape 
factor  using  marker  line  such  as  at  arrow. 
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Microcrack  «1  for  bOO  cycle  increments  in  fatigue  cycles 
from  2.5x103  to  b.OxlO^  cycles.  Crack  length  at  (a)  is 
30pm,  at  (b)  is  88pm. 
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fatigue  cycle  increments.  The  upper  tip  of  crack  #1  is  referred  to  as  Tip  A 
and  the  lower  tip  as  Tip  B.  The  average  rate  of  propagation  of  each  tip  of  #1 
with  fatigue  is  given  in  Fig.  21.  Figure  22  shows,  as  a function  of  fatigue, 
the  estimated  closure  stress  for  tips  A and  B of  crack  #1  and  Fig.  23  the 
crack  shape  factor.  The  occ  values  and  a/2c  were  obtained  using  the  proce- 
dures previously  described.  Averages  of  the  values  of  the  parameters  were 
calculated  between  two  successive  fatigue  increments  to  generate  the  data 
plotted. 

The  utility  of  a Paris  equation  corrected  for  crack  closure  [R12]  to 
describe  microcrack  propagation  is  accessed. 


W ■ A tSKeff]m  , (12) 

where 

aKeff  = (omax-occ)v/r/Q  F( a/c , <p ) . (13) 

A and  m are  material  constants.  Q is  a correction  for  plastic  zone  size, 
which  is  itself  dependent  upon  the  crack  shape  factor.  Using  the  notation  of 
Smith  and  Sorenson, [R20]  the  parameter  F(a/c,$)  is  used  to  describe  the 
variation  in  stress  intensity  range  along  the  crack  front.  The  case  of 
interest  here  is  $ = 90°,  corresponds  to  propagation  along  the  surface. 

Values  for  F(a/c,  41  = 90°)  vs  a/c  were  obtained  for  use  in  the  present 
calculation  by  extrapolating  to  shallow  crack  depth  values  for  F obtained  by 
Smith  and  Sorenson.  Thus,  the  approximation  has  been  made  that  the  cracks  are 
semi-elliptical,  planar,  normal  to  the  applied  stress  and  are  unaffected  by 
the  presence  of  intermetal  11c  particles  such  as  near  tip  B of  crack  #1. 

Predicted  to  measured  dc/dN  values  for  crack  #1  are  compared  in 
Fig.  24.  Here,  a value  of  m = 2 has  been  assumed,  typical  of  that  encountered 
for  A1  alloys  for  macrocrack  propagation  corrected  for  closure.  The  dotted 
line  is  the  value  of  dc/dN  which  would  correspond  to  a occ  = 0 case.  In 
Fig.  25,  dc/dN  is  compared  to  &Keff  for  the  four  test  cracks  studied,  from 
which  a best  linear  regression  fit  of  m * 1.6  is  obtained.  Figure  25  is  taken 
from  additional  analysis  of  the  microcrack  propagation  data  supported  by 
Rockwell  IR&D  funds. 
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Fig.  23  Crack  shape  factor  estimated  from  crack  compliance. 
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FATIGUE  CYCLES 

Fig.  24  Comparison  of  predicted  to  measured  cracking  rate  for  crack  #1. 
Dots  are  the  predicted  values.  Uppermost  curve  is  predicted 
cracking  rate  in  the  absence  of  crack  closure. 
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3.2.4  DISCUSSION 

From  experience  with  macrocrack  propagation  in  fracture  mechanics 
specimens,  one  generally  associates  a large  cyclic  range  In  CTOD  during 
loading  with  a correspondingly  large  range  in  stress  Intensity  and  conse- 
quently with  large  crack  propagation  ranges. [R7, 18]  This  model  requires 
modification  for  microcracks  because  the  plastic  zone  size  is  not  determined 
by  continuum  mechanics,  but  instead  is,  approximately,  given  by  the  size  of 
the  grain  at  the  crack  tip.[D8]  With  a large  grain  at  the  tip,  a large  CTOD 
can  be  achieved  with  a small  applied  stress  intensity.  The  accumulated  damage 
is  smaller  per  cycle,  compared  to  an  analogous  large  CTOD  in  a fracture 
mechanics  specimen,  because  the  strain  from  which  the  CTOD  develops  is 
distributed  over  a proportionally  larger  volume. 

To  appreciate  the  success  of  the  Akeff  analysis,  it  must  be 
recognized  that  there  is  no  meaningful  correlation  between  dc/dN  and  aK 
(uncorrected  for  closure).  Present  data  are  taken  very  near  the  threshold 
stress  intensity  range  value  for  propagation,  with  aK  < 2 MN/m3/2.  Crack 
closure  stress  is  therefore  a major  factor  in  determining  the  rate  of 
propagation  of  surface  microcracks  in  A1  2219-T851.  The  simple  Akef^  analysis 
also  provides  an  empirical  means  to  predict  the  cracking  rate  for  individual 
microcracks,  because  the  closure  stress  can  be  estimated  from  the  size  of  the 
grain  at  a crack  tip.[D8] 

Although  it  is  useful  in  this  case,  the  AKeff  description  of 
microcrack  propagation  is  a simplification  on  reality.  It  is  almost  certain 
that  the  cessation  of  reverse  yielding  in  the  plastic  zone,  on  unloading,  is  a 
function  of  location  in  the  zone  for  a given  external  stress  - not  the  sharp 
transition  which  the  present  model  assumes.  Further  analytic  insight  is 
necessary  for  special  problems  applicable  to  regime  where  microstructure  con- 
strains deformation.  It  is  probable  that  computer  modeling  of  closure  pro- 
cesses for  microcracks  would  be  highly  useful  in  development  of  an  improved 
effective  stress  intensity  range  description  of  the  early  propagation  of 
microcracks. 
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3.2.5  SUMMARY 

Microcrack  propagation  rates  are  measured  using  scanning  electron 
microscopy  for  surface  microcracks  developed  in  A1  2219-T851  by  fully  reversed 
loading  cyclic  fatigue.  Propagation  Is  followed  for  selected  cracks  In  the 
range  of  25  to  lOOym  In  surface  length.  At  Increments  in  the  fatigue  life- 
time, the  closure  stress  at  Individual  microcrack  tips,  at  the  surface,  is 
estimated  from  a measurement  of  the  crack  opening  displacement  at  the  tip. 
Crack  shape  factor  Is  estimated  from  a measurement  of  microcrack  compliance 
for  the  same  fatigue  Increments.  The  closure  and  crack  shape  data  are  used  to 
calculate  an  effective  stress  Intensity  range,  AKeff,  for  comparison  to 
measured  microcrack  propagation  rates.  The  results  Indicates  that  closure 
stress  is  a major  factor  In  determining  the  rate  of  propagation  of  surface 
microcracks  In  A1  2219-T851. 

3.3  THE  RELATIONSHIP  OF  CRACK  TIP  OPENING  DISPLACEMENT  TO  GRAIN 

SIZE  FOR  MICROCRACKS  IN  A1  2219-T851  [010] 

Departures  between  theoretical  and  experimental  values  of  plastic 
zone  size  have  recently  been  reported  for  low  values  of  stress  Intensity  in 
fracture  mechanics  specimens. [R18]  The  differences  are  attributed  to  a 
failure  of  the  continuum  theory.  A related  effect  has  also  been  observed  for 
surface  microcracks  with  lengths  of  approximately  the  grain  size  in 
A1  2219-T851.[D7,8]  In  this  instance,  analysis  of  microcrack  tip  opening 
displacement  (CTOD)  measurements  led  to  the  suggestion  that  the  zone  of 
significant  plastic  deformation,  for  the  special  case  of  microscopic  cracks, 
is  determined  by  the  size  of  a grain  at  the  surface  crack  tip.  This 
observation  Is  of  particular  interest  because,  potentially.  It  leads  to  a 
method  to  predict  the  effect  of  grain  size  on  the  rate  of  growth  of 
microcracks  during  fatigue.  It  has  been  previously  shown  that  the  magnitude 
of  crack  closure  resulting  from  residual  tensile  deformation  at  a microcrack 
tip  can  be  determined  by  measuring  CT0D.[D8]  Furthermore,  an  effective  stress 
intensity  range,  aK6^,  calculated  by  using  such  closure  values,  has  led  to 
adequate  prediction  of  rates  of  microcrack  growth  in  A1  2219-T851.[D9] 
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In  the  experiments  described  here,  correlations  were  sought  between 
CTOD  and  the  distance  of  the  crack  tip  to  the  grain  boundary  for  surface 
microcracks  In  A1  2219-T851.  Flexural  fatigue  specimens  of  an  alloy  with 
composition  6.3  wt%  Cu,  0.3%  Fe,  0.2%  Mn,  balance  aluminum  were  prepared  for 
minimum  residual  surface  stress  by  machining  with  sequentially  reduced  cutting 
depths.  The  specimens  were  polished  with  a 0.05pm  AI2O3  powder  and  then  given 
a light  chemical  etch  to  permit  Identification  of  the  grain  boundaries.  A 
triangular  shaped  constant  surface  stress  specimen  design  was  used,  mini- 
aturized to  permit  Insertion  Into  a scanning  electron  microscope  ( SEM) . [D4] 
Loading  of  a specimen  to  a desired  tensile  surface  stress  In  the  microscope 
was  accomplished  using  a jig,  described  elsewhere, [04]  which  displaced  the  end 
of  the  specimen  a measured  amount.  Specimens  were  fatigued  In  air,  for  either 
5 or  60%  relative  humidity  (RH),  with  the  principal  stress  axis  In  the  rolling 
direction  (for  which  <jy^d  = 350  MN/m2),  and  then  transferred  to  the  micro- 
scope for  study.  During  fatigue,  the  loading  was  fully  reversed  with  a peak 
cyclic  stress,  o,,^,  of  0.9  Oyje-|d.  Approximately  2500  fatigue  cycles  were 
applied  to  each  specimen  to  achieve  a distribution  of  microcrack  sizes  on  the 
surface. 

Figures  7 and  8 are  used  to  define  the  crack  measurement  nomenclature 
employed  here.  The  crack  length,  2c,  Is  measured  at  the  surface  from  tip  to 
tip.  The  distance  of  a crack  tip  to  the  next  grain  boundary,  cal'ed  z0.  Is 
measured  along  the  z axis,  which  Is  normal  to  the  stress  axis.  The  origin  of 
the  z axis  Is  placed  at  the  center  of  the  crack  at  the  surface,  a Is  the  sur- 
face stress  and  the  CTOD  measurements  were  made  for  0 ~ omax.  To  reduce  ad- 
ditional complexity  which  might  arise  In  Interpretation  of  how  to  measure  zQ 
with  variable  grain  shape,  measurements  were  restricted  to  grains  elongated  In 
the  rolling  direction,  with  crack  propagation  normal  to  the  rolling  direc- 
tion. Crack  opening  displacement  as  a function  of  location  along  a crack  Is 
represented  by  5(0, a)  and  was  measured  parallel  to  the  applied  stress.  Two 
crack  opening  displacement  values  were  found  for  each  crack  tip  studied; 

1)  An  opening  displacement  at  the  crack  center  termed  5(0, o);  2)  and  at  the 
crack  tip  termed,  6(c,o).  5 ( c ,0)  was  obtained  by  extrapolation  of  5 (z,o)  near 
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the  tip,  to  z = c.  6(c,o)/6(0,o)  has  previously  been  related  to  the  closure 
stress,  occ,  at  a microcrack  tip. 

The  ratios  5(c,o)/6(0,o)  and  z0/2c  are  compared  for  the  two  test 
humidities  In  Figs.  26  and  27.  z0  Is  normalized  to  the  microcrack  length 
because  It  Is  known  that  the  crack  closure  stress  Is  Independent  of  crack 
length  for  a given  surface  trace  of  the  cracking  path.[D4]  Each  data  point  Is 
from  an  Individual  microcrack  tip  of  a crack  for  which  30<2c<150um. 
6(c,o)/6(0,o)  Is  Interpreted  as  a measure  of  tensile  strain  at  the  crack 
tip.  The  minimum  value  of  6(c,o)/5(0,o)  ~ 0.1  and  corresponds  to  the  expected 
elastic  limit  of  CT0D.[R16]  Thus,  It  Is  found  that  6{c,o)/6(0,a)  Increases 
linearly  with  z0/2c  with  approximately  no  plastic  deformation  for  z0  * 0.  A 
shallower  slope  for  the  60%  RH  data  Is  attributed  to  a reduced  ductility  of 
the  material  at  the  crack  tip  as  a result  of  alteration  by  humidity.  Indeed, 
the  slope  of  the  6(c,o)/6(0,o)  vs.  z0/2c  line  can  be  Interpreted  as  a 
measurement  of  the  local  ductility  of  the  aluminum  matrix  In  the  microscopic 
region  at  the  crack  tip. 

Also  of  Interest  Is  to  use  the  data  of  Figs.  26  and  27  to  relate  the 
closure  stress  at  the  crack  tip  to  z0/2c;  where  crack  closure  Is  used  In  a gen- 
eric sense  to  refer  to  all  processes  leading  to  a reduced  cyclic  plastic  de- 
formation per  cycle  at  the  crack  tip,  and  hence  to  aKeff  « (omax"0cc^*  the 
maximum  cyclic  stress  minus  the  closure  stress.  For  surface  microcracks  In 
aluminum  alloys,  occ  can  arise  from  two  Independent  mechanisms:  (a)  physical 
contact  of  opposing  fracture  surfaces;  (b)  as  a result  of  residual  tensile 
strain  at  the  crack  tip.  Unlike  the  mechanism  (a),  often  observed  on  macro- 
cracks In  fracture  mechanics  specimens,  mechanism  (b)  Is  not  typically  found 
to  lead  to  physical  contact  behind  the  crack  tip  at  occ.  Instead,  “closure" 
at  a positive  external  load  on  unloading  Is  Interpreted  as  resulting  from  a 
cessation  of  reverse  yielding  In  the  plastic  zone.  Mechanism  (a)  is  of 
special  importance  In  the  very  early  stage  of  microcrack  growth,  just  as  a 
crack  leaves  an  Intermetalllc  particle  from  which  it  has  nucleated,  In  which 
case,  occ  can  approach  For  grain  sized  cracks,  surface  roughness  also 

places  a lower  limit  on  the  closure  stress  of  occ  ~ 0.2  for  A1  2219-T85I.  It 
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Fig.  27  Relationship  of  CTOD  to  distance  of  crack  tip  to  grain 
boundary  for  60%  RH. 
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is  the  mechanism  giving  rise  to  the  largest  occ  at  a point  on  a microcrack 
front  which  determines  occ  at  that  point. 

Prior  comparison  of  S(c,o)/S(0,o)  to  occ  [Fig.  14]  (Instead  of  z0/2c) 
has  found  a trend  similar  to  that  reported  In  Figs.  26  and  27,  differing  only 
In  behavior  for  small  6(c,o)/6(0,o).  CTOD  measures  only  the  effect  of  tensile 
residual  at  the  crack  tip  on  closure  and  If  <5(c,o)/6(0,o)  Is  small,  occ  Is 
determined  by  surface  roughness.  To  obtain  a relationship  between  occ  and 
z0/2c,  we  Ignore  data  for  occ  ~ 0.2  In  the  6(c,o)/5(0,o)  vs  occ  measurement. 
This  leads  to 


for  5%  RH,  and 


for  601  RH.  The  relationships  describe  only  the  effect  of  residual  tensile 
strain  at  the  crack  tip  on  the  closure  stress,  and  larger  values  of  occ  from 
roughness  may  supercede  that  given  by  the  relationships  in  special  circum- 
stances. Additionally,  the  data  from  Figs.  26  and  27  were  obtained  for  cracks 
which  had  propagated  beyond  the  grain  of  nucleatlon.  Within  the  grain  of  nuc- 
leatlon,  crystallographic  cracking  tends  to  decrease  the  closure  stress,  and 
It  Is  found  that  (14)  and  (15)  define  upper  bounds  of  the  observed  closure 
stress  In  that  case.  For  sufficiently  long  crack  length,  one  also  expects  the 
plastic  zone  size  to  revert  to  that  given  by  a continuum  analysis  for  the 
material  studied.  The  crack  size  at  which  this  occurs  Is  presently  known  only 
approximately  ( «•  400  um),  corresponding  to  an  average  of  four  grain  diameters 
In  the  alloy  Investigated.  As  a final  comment.  It  Is  of  particular  Interest 
that  measurement  of  zQ  at  the  surface  suffices  to  determine  trends  in  closure 
stress.  Undoubtly,  scatter  In  the  data  In  Figs.  26  and  27  relate  In  part  to 
the  crystallographic  orientation  of  Individual  grains  and  to  subsurface  grain 
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the  surface  Is  easier 
and  hence  closure 
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4.0  FATIGUE  BEHAVIOR  PREDICTIONS 

Models  of  nucleatlon  and  early  crack  growth  have  been  used  to  develop 
a computer  (Monte  Carlo)  simulation  of  fatigue  failure.  While  the  simulation 
Is  restricted  to  failure  arising  from  crack  nucleatlon  at  surface  particles  we 
believe  the  techniques  discussed  can.  In  the  future,  be  generalized  to 
describe  much  more  complex  failure  processes.  Crack  coalescence  remains  to  be 
the  least  understood  of  the  three  sequential  elements  of  failure  and  will 
require  substantial  amounts  of  further  investigation.  Presently,  coalescence 
enters  the  simulation  only  from  the  statistics  of  location  of  the 
microcracks.  We  have  observed  that  for  cracks  of  no  more  than  a few  grain 
diameters  In  length,  the  surface  tips  of  neighboring  cracks  do  not  sense  each 
others'  presence  until  the  tips  are  In  a common  grain.  This  can  be  expected 
from  the  nature  of  the  deformation  described  for  single  microcrack  growth. 

Once  in  a common  grain,  however,  changes  In  crack  tip  closure  stress  are  com- 
plex and  remain  to  be  described.  We  make  some  simplifying  assumptions  about 
the  nature  of  the  coalescence  process  once  crack  tips  are  in  a common  grain. 

Because  of  the  success  of  the  ONR  supported  research,  Rockwell  Inter- 
national is  continuing  support  of  the  lifetime  modeling  along  the  lines  first 
identified  in  the  ONR  effort.  In  addition  to  the  modeling  experiments  the  ONR 
funding  was  used  to  develop  the  one-dimensional  simulation  of  fatigue  failure 
described  In  Section  4.1.  The  experimental  data  against  which  the  predictions 
are  tested  were  obtained  with  Rockwell  support.  Subsequent  to  the  completion 
of  the  ONR  effort,  a more  general  two-dimensional  model  of  nucleation  has  been 
developed  and  we  discuss  this  briefly  to  illustrate  potentials  for  research  in 
this  area  in  the  future.  Section  4.1  is  comprised  of  comparison  of  predicted 
to  measured  fatigue  behavior  on  a microscopic  level  and  therein  is  described 
the  general  simulation  approach.  In  Section  4.2  further  details  of  the  crack 
coalescence  simulation  are  provided  along  with  predictions  of  alloy  and 
humidity  effects  on  scatter  and  mean  values  of  fatigue  lifetime. 
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4.1  A MONTE  CARLO  SIMULATION  OF  MICROCRACK  GROWTH 

A Monte  Carlo  based  computer  simulation  of  the  fatigue  crack  Initia- 
tion phase  of  an  aluminum  alloy  is  described.  Processes  modeled  Include 

1)  nucleatlon  of  microcracks  at  Intermetallic  sites  at  the  alloy  surface;  and 

2)  early  microcrack  propagation  wherein  cracking  rate  is  retarded,  and  highly 
sensitive  to  the  alloy  microstructure  In  the  region  of  a microcrack.  The 
simulation  is  used  to  predict  the  distribution  In  closure  stress  values  for 
surface  microcracks  as  a function  of  crack  length  and  the  number  of  fatigue 
cycles.  Predicted  and  measured  distributions  In  closure  stress  values  are 
compared  using  A1  7075-T7  alloys  for  two  different  grain  sizes. 

4.1.1  INTRODUCTION 

Analytical  models  of  nucleatlon  and  early  growth  of  microcracks  on 
the  surface  of  aluminum  alloys  subjected  to  fatigue  loading  have  recently  been 
proposed  and  evaluated. [D2, 8-10]  The  models  are  simplistic  and  yet  contain 
sufficient  reality  to  make  reasonable  predictions  of  the  rates  of  crack 
formation  and  propagation  in  the  initiation  stage  of  the  fatigue  failure 
process.  We  describe  here  an  extension  of  this  research,  in  which  the 
microscopic  failure  process  models  are  used  to  develop  a Monte  Carlo 
simulation  of  fatigue  failure  of  certain  aluminum  alloys.  The  general 
approach  is  to  nucleate  cracks  in  simulation,  which  then  propagate  according 
to  rules  relating  cracking  rate  to  the  location  of  the  surface  microcracks 
relative  to  neighboring  grains.  In  more  advanced  versions  of  the  simulation, 
crack  coalescence  is  allowed.  Specimen  fatigue  lifetime  is  determined  by  the 
smallest  number  of  fatigue  cycles  required  for  any  crack  present  on  a specimen 
to  reach  a critical  size.  Scatter  In  fatigue  lifetime  from  (simulated) 
specimen  to  specimen  then  results  from  statistical  fluctuations  In  alloy 
microstructure  consistent  with  the  average  alloy  microstructure. 

An  alternate  form  of  such  a simulation  Is  to  predict  distributions  in 
microcrack  parameters  with  fatigue.  This  can  Include  numbers  of  cracks,  dis- 
tribution in  crack  size,  and  parameters  of  Importance  to  remaining  lifetime 
prediction  - such  as  distribution  In  crack  closure  stress  with  microcrack 
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length.  Such  microscopic  predictions  have  several  applications  Including:  (a) 
definition  of  fatigue  damage  at  Intermediate  stages  of  failure  - for  use  In 
analyzing  nondestructive  evaluation  techniques  to  determine  remaining  fatigue 
lifetime,  (b)  to  provide  simulated  microcrack  parameters  for  comparison  to 
measured  parameters.  It  Is  this  latter  use  that  Is  of  interest  here,  as  It 
provides  a means  to  determine  the  strengths  and  weaknesses  of  the  models  on 
which  the  simulation  Is  based. 

Specifically,  we  discuss  a comparison  between  predicted  and  measured 
values  of  the  average  closure  stress  of  surface  microcracks,  as  a function  of 
crack  length,  for  two  A1  7075  17  alloys.  Crack  closure  is  a prime  factor  In 
determining  the  rate  of  propagation  of  the  surface  microcracks  during 
fatigue.  An  Important  mechanism  of  closure  stress  development  is  from 
residual  tensile  strain  at  the  microcrack  t1p.[D8]  It  has  previously  been 
shown  that  when  this  mechanism  pertains,  the  magnitude  of  the  closure  stress 
is  sensitive  to  the  size  of  the  grains  ahead  of  the  propagating  micro- 
crack. [010]  Consequently,  the  simulated  versus  experimental  comparisons  are 
made  for  two  alloys  having  different  mean  grain  sizes.  In  the  direction  of 
crack  propagation  (75  urn  and  25  pm). 

While  we  do  not  report  on  predicted  fatigue  lifetime  In  the  present 
section  the  general  form  of  the  simulation,  used  with  respect  to  lifetime 
prediction.  Is  described,  as  the  microscopic  predictions  discussed  were 
extracted  as  a by-product  of  the  simulation.  The  simulation  technique  de- 
scribed Is  an  early  version  of  a program  which  has  subsequently  been  revised 
to  reflect  further  mechanistic  Insights  and  to  Incorporate  additional 
sophistication  In  approximations  required  In  the  representation  of  reality. 
Limitations  and  approximations  In  the  simulation,  used  to  generate  the 
predictions,  are  discussed  as  well  as  opportunities  and  approaches  for 
Improvement  of  the  simulation. 
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4.1.2  THE  FATIGUE  MODEL 

We  begin  by  describing  the  sequence  of  fatigue  failure  events  to  be 
modeled  and  the  general  approach  used  to  simulate  fatigue  failure,  and  then 
discuss  the  specific  failure  process  models  employed  and  the  approximation 
made  encoding  them  Into  a computer  simulation  of  microcrack  growth. 

Microscopic  Failure  Sequence 

Three  major  processes  can  Influence  the  fatigue  lifetime  and  scatter 
In  lifetime  of  structural  alloys:  (1)  microcrack  nucleatlon;  (2)  early  micro- 
crack growth,  and  (3)  microcrack  coalescence.  These  sequential  events 
comprise  the  Inltatlon  stage  of  the  fatigue  failure  process  and  precede  the 
formation  of  a terminal  macrocrack  whose  subsequent al  progress  can  be 
predicted  by  fracture  mechanics.  Our  modeling  Is  presently  constructed  to 
describe  fatigue  failure  wherein  crack  nucleatlon  Is  at  or  near  the  surface  at 
particles.  It  Is  therefore,  an  analog  of  failure  of  peak  hardness  or  overaged 
2000  and  7000  series  aluminum  alloys,  which  contain  Intermetalllc  particles  of 
sufficient  size  to  enhance  nucleatlon  rates  substantially  above 
crystallographic  or  at  grain  boundary  cracking.  The  simulation  describes  a 
material  relieved  of  long  range  residual  surface  stress,  having  no  major 
subsurface  defects,  and  currently  also  Ignores  the  effect  of  trlaxlally  of  the 
applied  stress  on  the  Initiation  process.  With  no  gross  subsurface  defects 
and  no  macroscopic  residual  surface  stress,  the  absence  of  a bulk  constraint 
of  microplastic  deformation  during  fatigue  ensures  that  crack  nucleatlon 
occurs  at  or  near  the  surface. 

Fatigue  failure  Is  a statistically  Influenced  process  and  the  scatter 
In  lifetime  can  arise  from  fluctuations  In  each  of  the  three  microscopic  types 
of  failure  events  comprising  initiation.  One  expects  the  relative 
contribution  of  the  three  processes  to  scatter  In  lifetime  to  be  principally  a 
function  of  the  number  and  density  of  microcracks  on  a specimen.  Thus  if  the 
cyclic  stress  amplitude  Is  sufficiently  small  (high  cycle  fatigue)  few  cracks 
are  formed,  nucleatlon  will  dominate  the  statistics  of  failure.  For  low  cycle 
fatigue,  with  multiple  crack  coalescence  commonplace,  the  statistics  of 
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coalescence  is  of  prime  Importance.  At  an  Intermediate  stress  statistical 
scatter  In  microcrack  growth  can  be  expected  to  govern  the  statistics  of 
failure. 

General  Simulation  Technique 

The  simulation  deals  only  with  the  growth  of  microcracks  as  seen  from 
the  surface.  A fixed  crack  shape  factor  (depth/length)  of  0.35  Is  assumed  for 
all  microcracks.  We  begin  by  defining  a function  P(N),  which  Is  the 
probability  of  nucleation  of  Isolated  cracks  at  Intermetalllc  particles  In  an 
element  of  surface  area,  S,  during  the  Nth  fatigue  cycle.  P( N)  Is  calculated 
for  a specific  cyclic  stress  and  specific  orientation  of  the  material  relative 
to  the  surface  and  to  the  principal  stress  axis.  For  the  case  of  flexural 
fatigue,  used  In  the  present  work,  the  stress  and  material  orientations  are 
constant  over  the  entire  gauge  section  of  the  specimen,  and  so  S Is  simply  the 
surface  area  of  the  specimen.  Calculation  of  P ( N)  for  more  complex  components 
requires  that  the  surface  be  divided  into  elements  of  approximately  constant 
material  orientation  and  cyclic  stress  amplitude. 

For  each  specimen  fatigue  lifetime  to  be  simulated,  P(N)  Is  used  with 
a random  number  selection  procedure  to  define  the  cycles  N^  at  which  the  1th 
crack  nucleates  with  a sequential  ordering  so  that  N.j+1  > N.j . Next,  for  each 
nucleation  event,  we  determine,  by  random  selection,  the  width  of  the 
Intermetalllc  from  which  nucleation  takes  place  (Fig.  3),  the  width  of  the 
grain  containing  the  Intermetalllc,  G^,  and  the  distance  of  the  intermetalllc 
from  the  grain  boundary,  Z^.(F1g.  28)  There  are  specific  constraints  on  the 
range  of  values  of  W^,  G^,  and  which  can  occur  for  a given  N which  arise 
from  the  mechanism  of  nucleation  as  discussed  later.  Change  In  microcrack 
length  with  Increment  In  fatigue  cycles  aN  Is  described  by 

AC  = AN  A (AKeff)m  , (16) 

where  c Is  the  half  crack  length,  m and  A are  material  constants  and  AKeff  is 
an  effective  stress  intensity  range.  AKe^f  <*  (omax  - occ),  the  maximum  cyclic 
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Fig.  28  Schematic  Illustration  of  surface  section,  defining  micro- 
structural  parameters  important  to  crack  nucleation  and 
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stress  minus  the  crack  closure  stress.  occ  is  determined  from  rules  which 
relate  the  closure  stress  to  the  Instantaneous  location  of  a microcrack 
relative  to  surrounding  microstructure. 

Simulation  of  growth  of  the  1th  crack  begins  for  N = with  the 

crack  length  2c  * W^.  An  Incremental  change  In  2c  is  calculated  for  each  half 
of  the  crack  using  Eqn.  (16).  The  new  2c  Is  used  to  determine  a new 
for  N * N + aN,  and  the  simulation  of  propagation  continues  In  that  manner. 

For  crack  1*1,  growth  Is  continued  until  2c  reaches  a critical  size  cr  at 
which  point  we  set  the  cycles  to  specimen  failure  Np  = N.  The  second  crack 
then  Is  followed  beginning  at  N2  cycles  and  N Is  compared  to  Np,  and  2c  to 
cr.  If  2c  reaches  cr  with  N < Np  the  new  Np  becomes  the  N value  of  the  second 
crack.  The  above  procedure  Is  terminated  when,  for  the  1th  crack,  N^  < Np. 
During  crack  growth  the  size  of  grains  neighboring  the  nucleatlon  site  are 
determined  as  needed,  just  as  propagation  of  a crack  tip  Into  a new  grain 
begins.  The  width  of  the  new  grain  is  selected  at  random  using  the 
distribution  In  width  of  grains  normal  to  the  stress  axis.  Crack  coalescence 
is  handled  In  much  the  same  manner,  by  first  defining  the  probability  function 
of  nucleatlon  of  cracks  In  grains  of  close  proximity.  The  cycles  to 
nucleatlon  of  each  crack  Is  selected  along  with  kthe  initial  parameters 
describing  the  nucleatlon  site  sizes  and  locations.  Propagation  of  each  crack 
Is  then  simulated  to  the  point  of  coalescence,  and  beyond  as  necessary. 

Single  Crack  Nucleatlon  and  Growth 

In  the  present  paper  we  compare  predicted  to  measured  parameters  for 
Isolated  cracks  and  hence,  defer  discussion  of  simulation  of  crack  coalescence 
to  later  reports.  The  nucleatlon  phenomenon  is  best  simulated  with  a two- 
dimensional  model.  For  this  initial  effort,  however,  we  chose  to  use  a one 
dimensional  simulation  of  the  nucleatlon  process.  This  provided  substantial 
savings  in  computation  time.  We  rationalized  the  approximation  on  the  basis 
that  the  closure  stress  distribution  of  Interest  Is  especially  sensitive  to 
the  crack  growth  rather  than  to  the  nucleatlon  simulation.  For  purposes  of 
comparison,  we  describe  the  approaches  for  both  the  one-  and  two-dimensional 


69 


Rockwell  International 

Science  Center 


SC5050.1FR 


simulations  of  crack  nucleation,  which  aids  in  definitions  of  simplifying 
approximations  made  In  the  one-dimensional  simulation.  We  begin  with  a 
general  description  of  the  nucleation  process  to  be  modeled. 

Crack  Nucleation:  We  use  as  a basic  equation  describing  cycles  to 
nucleation  at  a surface  particle  an  expression  proposed  by  Chang  et.al.[R2] 


^max  * V 

With  reference  to  Fig. 28,  N is  the  number  of  cycles  to  onset  of  crack 
propagation  from  a particle  into  the  matrix,  W Is  the  width  of  the  particle 
normal  to  the  stress  axis,  D is  the  largest  surface  slip  distance  (at  45°  to 
the  principal  stress)  within  the  grain  containing  the  particle,  omax  is  the 
peak  cyclic  stress  amplitude,  and  C0  and  a0  are  material  constants.  The 
expression  pertains  to  fully  reversed  constant  amplitude  loading,  and  with 
simplifying  assumptions  describes  the  effect  of  increasing  strain  energy 
density  at  a nucleation  site  which  accrues  with  cyclic  microplastic 
deformation.  Equation  (17)  defines  a lower  bound.  In  that  for  a particle  with 
a given  W,D  and  omax  it  predicts  the  minimum  value  of  N at  which  cracking  into 
the  matrix  will  occur.  Crystallographic  orientation  of  the  grains  was  ignored 
in  derivation  of  (17)  and  less  favorable  orientations  are  one  factor  leading 
to  scatter  of  nucleation  to  large  N values. 

Consider  the  progression  of  nucleation  one  typically  observes  with 
A1  2219-T851,[D12]  wherein  the  nucleation  sites  are  at  intermetalllc 
particles.  During  fatigue  the  first  nucleatlons  occur  at  the  largest 
intermetal  lies  located  in  large  grains  permitting  the  largest  possible  values 
of  D ( Dmax ) * The  possible  cycles  to  the  first  nucleation  is  Nmin  * 

C0/(  J Wmax  Dmax  - o0)2).  As  N increases,  additional  nucleation  events 

occur  for  Intermetal  1 ics  corresponding  to  progressively  smaller  values  of  W 
and  D.  For  a given  N the  smallest  intermetalllc  permitting  nucleation  is  W^p 
* Cq2/(N  Dmax(amax  - o0)2)2.  One  actually  observes  nucleation  for  all 
particle  sizes  Wmax  > W > Wmin.  This  is  in  part  because  crystallographic 
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orientation  effects  nucleatlon  at  surface  particles,  but  is  also  because 
variation  In  the  internal  structure  from  one  Intermetal  lie  to  the  next 
Influences  nucleatlon.  The  percentage  of  specimen  surface  area  in  which 
nucleation  is  possible  Is  Identically  zero  at  N = N^,,,  and  Increases  with 
Increasing  N - Initially  as  portions  of  the  largest  grains  In  the  alloy 
neighboring  the  grain  boundaries  become  potential  areas  for  nucleatlon. 


Two-Dimensional  Nucleatlon  Expression:  The  crack  nucleatlon 
phenomena  we  wish  to  simulate  Is  most  accurately  described  by  a two- 
dimensional  model.  For  point  of  reference,  we  describe  briefly  a two- 
dimensional  procedure  to  determine  the  number  of  fatigue  cycles  for  successive 
nucleation  of  cracks  on  the  surface  of  a simulated  test  specimen.  We  assume  a 
random  distribution  of  Intermetallic  sizes  over  the  surface.  The  probability 
of  nucleatlon  per  fatigue  cycle  P( N)  is,  then,  calculated  from  the  distribu- 
tion in  possible  slip  distances,  A0(D),  and  the  distribution  In  Intermetallic 
sizes,  F(W).  F(W)  Is  the  number  of  intermetal  lies  per  unit  area  with  width  W 
normal  to  the  principal  stress.  A+45°(D)  and  A_4go(D)  are  normalized  (to 
unity)  relative  probability  distributions  of  finding  a cord  length,  D,  across 
a grain  at  +45°  and  -45°  from  the  stress  axis  respectively.  P(N)  is  simply 
taken  to  be  proportional  to  the  number  of  Intermetal  lies  available  for 
fracture  as  defined  by  the  limits  of  (17). 


P(N)  =SC 


i 

/ 


max 


F(W) 


Gf+45,-45  Wmin(N) 


r 


MX  (D  - Dm(n) 


min(N.W) 


A0(D)  dD  dW. 


(18) 


In  writing  (18)  we  have  assumed  that  there  Is  no  substantial  saturation  of 
nucleatlon  for  any  Intermetallic  size,  i.e.,  that  few  particles  of  any  given 
size  fracture  during  fatigue  - even  those  that  are  located  in  grains  having 
the  most  favorable  orientations.  S Is  the  specimen  surface  area,  Cp  is  the 
probability  of  nucleatlon  per  cycle  at  a particle  meeting  the  condition  for 
nucleatlon  defined  by  (17).  Currently,  Cp  must  be  determined  experimentally. 

Dm1n  = co^N^^°max  ” °o^‘ 
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One-Dimensional  Nucleaton  Expression:  Using  the  same  notation  as 
above  we  set 


PIN) 


■scp  / 


max 


AID)  dD 


Jm1n(W  * 20  um,N) 


pm 

10  uin 


f (W)  dW 


(19) 


We  take  AID)  to  be  the  normalized  distribution  In  grain  widths  normal  to  the 
principal  stress  axis  - which  Is  equivalent  to  a simulation  In  which  the 
grains  are  equiaxial  In  the  plane  of  the  surface.  In  calculating  Dmirl  we  use 
a fixed  Intermetal  11c  width  of  W = 20  urn,  corresponding  to  an  average  of  the 
size  of  intermetalllc  for  which  nucleatlon  occurs  In  the  7075  alloys  to  be 
simulated.  Using  (19)  the  cycles  Nj,  at  which  a sequence  of  cracks  1 = 1,300 
are  nucleated,  are  selected  randomly. 

The  width  of  the  grain,  Gj,  In  which  the  1th  crack  nucleates  Is  then 
chosen  using  AID)  for  D > Dm1n  (w  = 20  ym,  N = Nj ).  We  define  a random  number 
0<:  R « 1, 


and  by  Inversion  a function 


G = X(R), 


(20) 


(21) 


from  which  G^s  chosen  by  random  selection  of  R.  Finally  the  location  of  the 
intermetalllc,  Z^,  within  the  grain  Is  chosen  at  random,  by  assuming  any 
location  In  the  grain  of  initiation  is  equally  likely. 
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With  respect  to  a two-dimensional  simulation  the  one-dimensional 
approximation  has  the  following  defects. 

1)  It  overestimates  the  absolute  numbers  of  nucleations  which 
occur,  although  not  necessarily  the  relative  change  in  P ( N ) 
with  N. 

2)  For  N = Nm^n  it  permits  nucleation  in  the  center  of  large 
grains  where  they  are  only  observed  at  the  periphery  of  the 
grains. 

3)  An  average  value  of  intermetal  lie  width,  and  hence,  initial 
crack  size  is  used. 

Early  Crack  Growth;  The  initial  crack  length  (2c)  is  set  to  20  pm 
for  all  nucleation  events.  For  the  ith  crack,  grain  width  G-j  and  location  of 
the  nucleation  site  Zn-  are  defined  as  described  above.  If  propagation  into  a 
neighboring  grain  occurs  during  growth  of  any  microcrack,  the  neighboring 
grain  width  is  selected  at  random  from  A(D),  (i.e.  the  width  of  the  grain  is 
found  using  (21)  for  Dm1n  = 0 in  (20). 

Change  in  crack  length  versus  aN  for  each  crack  is  found  using  (16) 
with  the  deterministic  parameter  being  the  closure  stress  occ,  used  to  calcu- 
late AKeff.  We  recognize  three  mechanisms  which  can  give  rise  to  acc,  of 
which  residual  tensile  strain  at  the  microcrack  tip  is  of  most  importance,  and 
is  the  only  mechanisms  simulated  here.[D8]  The  neglected  mechanisms  involve 
contact  of  opposing  fracture  surfaces,  and  pertain  only  when  acc  is  small  or 
if  2c  = W.  Residual  tensile  strain  at  the  crack  tip  gives  rise  to  acc  de- 
scribed by: 


cc 


a zc 
~ZT 


F(R) 


max 


for  cracks  within  the  grain  of  nucleation,  and 


(22) 
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_ “ Zo 

°cc  2c  amax 


(23) 


for  cracks  beyond  the  grain  of  nucleation.[D10]  a is  a material  constant.  z0 
is  the  distance  of  a crack  tip  to  the  next  grain  boundary  as  illustrated  by 
Fig.  7.  Function  F(R)  has  values  lying  between  zero  and  unity. 

Thus,  as  represented  by  Eqns.  (22)  and  (23)  the  crack  closure  stress 
varies  with  microcrack  propagation  with  largest  values  encountered  as  a crack 
tip  enters  into  a new  grain.  F(R)  is  included  to  account  for  the  fact  that 
crystallographic  propagation,  which  more  likely  occurs  within  the  grain  in 
which  a crack  nucleates,  is  often  associated  with  lower  values  of  closure 
stress.  For  the  7075  alloys  to  be  modeled  we  set  F(R)  = R where  0 < R < 1.  A 
final  constraint  used  on  Eqns.  (22)  and  (23)  occurs  if  the  calculated  value  of 
°cc  > °max*  case  we  set  OqC  - ^max* 

To  begin  the  growth  simulation  for  the  ith,  crack  F(R)  is  selected  at 
random  and  <jcc  is  calculated  for  each  of  the  two  surface  crack  tips  using 
Eqn.  (22).  In  the  event  that  the  intermetallic  overlays  a grain  boundary  F(R) 
is  found  independently  for  each  crack  tip.  Once  propagation  into  a 
neighboring  grain  begins,  Eqn.  (23)  is  used  to  calculate  acc.  The  values  of 
occ  are  updated  for  each  aN,  as  2c  changes  with  microcrack  growth. 


4.1.3  EXPERIMENTAL  PROCEDURES 

The  two  A1  7075-T7  alloys  utilized  were  taken  from  the  same  heat  of 
material  having  a composition  wt  % 1.6  Cu,  2.5  mg,  5.6  Zn,  0.3  Cr,  balance 
Al.  The  large  grained  alloy  was  in  a T731  condition.  The  small  grained  alloy 
was  prepared  by  a thermal  mechanical  process.  A(D)  distributions  in  the  long 
transverse  direction  (the  direction  of  crack  propagation)  for  the  two  alloys 
are  shown  in  Fig.  29.  The  material  constant,  a,  for  use  In  Eqns.  (22)  and 
(23)  was  obtained  by  comparing  crack  tip  opening  displacement  to  distance  of 
the  crack  tip  to  the  grain  boundary  using  the  large  grained  material .[DIO] 
Flexural  fatigue  specimens,  with  a tapered  design  to  achieve  constant  surface 
stress,  were  manufactured  from  the  center  of  a 1.0  cm  thick  plate  by  machining 
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with  decreasing  cutting  depths  to  minimize  residual  stress  at  the  surface. 
Speciments  were  then  mechanically  polished  using  a 0.05  um  AI2O3  powder. 

Fatigue  loading  was  fully  reversed  with  °max  = °*9  ayield»  in  60% 
relative  humidity  air  and  with  the  stress  axis  parallel  to  the  rolling 
direction;  It  was  terminated  when  maximum  surface  microcrack  lengths  of 
approximately  200  yin  were  achieved.  The  flexural  fatigue  geometry  employed  is 
designed  for  insertion  Into  a scanning  electron  microscope,  and  using  a jig  in 
the  SEM  the  specimen  can  be  loaded  to  a desired  tensile  surface  stress  by 
deflecting  the  specimen  end  a measured  amount. 

Using  a scanning  electron  microscope,  the  average  closure  stress  of 
selected  surface  cracks  was  measured  as  a function  of  microcrack  length.  The 
measurement  techniques  have  been  discussed  before.  Briefly,  they  rely  upon 
measurement  of  two  values  of  crack  opening  displacement,  from  which  an 
empirical  estimate  of  the  average  closure  stress,  acc,  is  made.[D4] 

4.1.4  RESULTS  AND  DISCUSSION 

The  computer  generated  predictions  of  occ  for  the  two  grain  sizes  are 
shown  in  Fig.  30  and  the  experimentally  measured  values  In  Fig.  31.  The  pre- 
dicted and  experimental  values  are  for  Isolated  cracks  developed  by  fatigue 
for  = 0.9  oyfeid*  occ  is  the  average  of  the  occ  values  of  the  two  sur- 
face crack  tips  for  each  simulated  crack.  A general  trend  towards  a smaller 
acc  for  the  fine  grained  material  Is  evident  and  is  predicted  by  the  model. 

It  Is  also  apparent  that  the  simulation  shows  Inadequacies  of  the  models  In 
three  areas. 

a)  For  short  crack  length  the  simulation  predicts  the  occur- 
rence of  more  cracks  with  large  values  of  6CC  than  are 
actually  observed.  We  believe  the  modeling  problem  is  that 
closure  stress  is  only  known  to  affect  Mode  I propagation. 
We  have  observed  that  shear  mode  propagation  continues  Into 
grains  for  which  large  closure  stresses  preclude  Mode  I 
propagation.  As  a consequence,  microcracks  are  not 
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Fig.  30  Predicted  values  of  occ  vs.  microcrack  lengths  for  two  grain 

sizes.  Each  data  point  is  from  a different  simulated  microcrack. 
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Fig.  31  Measured  values  of  occ  vs.  microcrack  length  for  two  grain  sizes 
Dash  lines  are  upper  bounds  of  the  data  from  the  simulation. 
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necessarily  permanently  trapped  In  a high  closure  stress 
condition  as  assumed  in  the  simulation. 

b)  We  expect  that,  for  sufficient  microcrack  length,  Eqn.  (23) 
will  fail  when  the  crack  length  becomes  so  large  that  a 
continuum  analysis  is  appropriate  for  pre  prediction  of  the 
plastic  zone  size.  The  occurrence  of  cracks  with  occ 
larger  than  expected  for  the  fine  grained  material  for  2c  > 
lOOum  Is  attributed  to  such  an  effect. 

c)  Experimentally,  we  observe  no  cracks  for  which  ?cc  < 0.15 

°max*  Presumed  t0  be  result  of  a contribution 

of  fracture  surface  roughness  to  closure,  which  sets  a 
minimum  value  to  acc,  and  which  has  been  neglected  in  the 
present  simulation. 

The  dashed  upper  bounds  of  the  occ  values  in  Fig.  31  are  the  upper 
bounds  from  the  simulation  and  suggest  that  the  general  trend  in  closure 
stress  with  crack  length  has  been  successfully  modeled.  Empirical  modifica- 
tion can  be  made  Immediately  to  refine  the  models  to  account  for  Items  (b)  and 
(c)  above.  With  respect  to  Item  (a)  further  research  Is  required  to  establish 
if  humidity  and  the  presence  of  a Mode  I closure  stress  affects  Mode  II 
microcrack  propagation. 

4.1.5  SUMMARY 

A Monte  Carlo  based  computer  simulation  of  fatigue  crack  initiation 
processes  pertinent  to  the  fatigue  failure  of  aluminum  alloys  is  discussed. 

The  simulation  Is  based  upon  models  which  describe:  (1)  nucleation  of  cracks 
at  intermetallic  particles  at  or  near  the  surface,  and  (2)  the  early  stage  of 
microgrowth.  An  effective  stress  intensity  range,  AKeff,  Is  used  to  calculate 
the  rates  of  propagation  of  individual  cracks.  In  turn,  the  microcrack 
closure  stress,  used  to  determine  Akeff  for  a given  microcrack,  is  determined 
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from  a model  predicting  the  effect  of  crack-tip  - to  grain  boundary  distance 
on  closure  from  residual  tensile  strain  at  the  crack  tip. 

The  simulation  is  applied  to  prediction  of  the  distribution  in  the 
average  closure  stress  of  microcracks  as  a function  of  microcrack  length  for 
two  A1  7075-T7  alloys;  a fine  grain  material  with  an  average  size  of  25  urn  and 
a large  grain  material  having  a 75  urn  average  grain  size  in  the  direction  of 
crack  propagation.  Simulation  predictions  are  compared  to  experimentally  de- 
termined closure  stress  distributions  for  the  two  alloys.  The  general  trends 
of  decreasing  closure  stress  with  microcrack  length  for  both  materials,  and  a 
smaller  average  closure  stress  for  the  fine  grained  material  is  predicted  by 
the  simulation.  The  experimental  data  also  indicate  three  areas  In  which  the 
simulation  needs  Improvement: 

1)  To  allow  for  Mode  II  crack  propagation  for  large  closure 
stresses  which  limit  Mode  I propagation  rates. 

2)  To  describe  the  transition  from  a plastic  zone  size 
determined  by  grain  size  for  smaller  microcracks,  to  one 
determined  by  a continuum  analysis  for  sufficiently  large 
cracks. 

3)  To  account  for  presence  of  a lower  bound  on  crack  closure 
stress  resulting  from  a mechanism  of  contacting  the 
opposing  fracture  surfaces  behind  the  crack  tip. 

4.2  ADDITIONAL  MODELING  CONSIDERATIONS 

Humi di ty : Relative  humidity  enters  the  simulation  both  in  the  nu- 
cleation  and  early  microcrack  propagation  processes.  One  major  effect  of 
humidity  is  to  decrease  the  ductility  of  a thin  layer  of  material  (=5  ym)  near 
the  surface.  This  occurs  as  a consequence  of  cyclic  loading  and  marks  a work 
hardening  which  accompanies  a penetration  of  some  as  yet  undefined  environ- 
mental species  into  the  surface,  presumably  by  a dislocation  sweep-in 
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mechanism.  The  reduced  surface  ductility  leads  to  a reduced  probability  of 
nucleation  reflected  by  changes  in  material  constants  C0  and  a0  in  the 
nucleation  Eq.  (17),  and  also  in  reduced  ductility  at  a crack  tip  and  hence 
smaller  values  of  a in  Eq.  (23).  These  two  effects  are  competing  and  one 
expects  the  net  effect  of  humidity  on  fatigue  lifetime  of  aluminum  alloys  will 
be  dependent  upon  the  distribution  of  Intermetallic  and  grain  sizes,  as  well 
as  of  the  direct  effect  of  humidity  on  the  material  coefficients. 

Several  additional  effects  of  humidity  must  be  considered  in  a com- 
prehensive model:  1)  Increased  humidity  Is  known  to  lead  to  an  unusual  weak- 
ness between  intermetal lie-matrix  1nterfaces,[Rl]  which  can  further  modify  the 
behavior  described  by  Eq.  (17);  2)  For  cracks  in  fracture  mechanics  specimens, 
the  environment  alters  the  closure  stress  only  near  the  crack  tip.  The  duc- 
tility of  the  material  at  the  crack  tip  changes  with  increased  humidity 
leading  to  smaller  occ.  However,  the  cyclic  deformation  deep  within  the 
plastic  zone  occurs  in  the  absence  of  environmental  modification  and  it 
appears,  at  least  for  certain  aluminum  alloys,  it  suffices  to  describe  the 
effect  of  humidity  using  Eq.  (16)  by  only  changing  acc,  with  material  con- 
stants A and  m Independent  of  humidity.  For  microcracks,  at  least  at  the 
surface,  the  entire  zone  of  deformation  is  affected  by  the  environment  and 
hence  one  can  expect  that  A,  m,  and  occ  will  all  be  a function  of  humidity, 
reflecting  not  just  the  effect  of  reduced  ductility  on  occ  but  also  a change 
in  vulnerability  of  the  material  to  cyclic  damage.  Presently,  the  various 
material  coefficients  and  their  altered  values  with  humidity  must  be  deter- 
mined experimentally  before  the  simulation  can  proceed.  If,  however,  in  the 
future  the  coefficients  are  related  to  alloy  composition  and  heat  treatment, 
the  fundamental  predictive  power  of  the  failure  simulation  will  be  substan- 
tially increased. 

Crack  Coalescence  In  Low  Cycle  Fatigue 

In  low  cycle  fatigue  crack  nucleation  density  is,  typically,  suffici- 
ent that  coalescence  of  microcracks  becomes  probable.  The  effect  of  coales- 
cence is  to  reduce  fatigue  lifetime  from  that  realized  with  Isolated  cracks. 
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The  importance  of  an  individual  coalescence  event  depends  upon  the  relative 
location  of  the  cracks  involved  including  the  sizes  of  the  grains  in  the 
region  containing  the  coalescing  cracks.  Presently  we  have  no  model  for  the 
changes  in  cyclic  stress  intensity  range  which  occurs  when  tips  of  microcracks 
enter  a common  grain.  From  fracture  mechanics  one  expects  an  increase  in 
cyclic  stress  intensity  range.  However,  we  know  that  counteracting  this 
effect  is  an  Increase  in  occ  associated  with  larger  tensile  strain  at  the 
crack  tips.  In  the  simulation  we  make  the  approximation  that  coalescing 
cracks  do  not  sense  each  other  until  their  surface  tips  touch,  at  which  point 
coalescence  occurs  (neglecting  the  effect  of  subsurface  crack  path).  Pres- 
ently we  deal  only  with  the  coalescence  of  two  cracks,  but  the  procedure  de- 
scribed below  can  be  extended  to  treat  multiple  cracks. 


Using  the  notation  of  Section  4.1  we  define  the  probability  of  nu- 
cleation  of  a second  crack  In  a grain  neighboring  a grain  containing  a crack 
developed  previously.  In  the  Nth  cycle,  over  the  surface  area  (S)  element,  to 
be 


P“(N) 


_ 2 P ( N ) 
' ~G5 


P(N)  dN 


(24) 


where  G is  the  number  of  grains/area.  The  total  number  of  single  cracks  in  S 
is  simply  N./~NP(N)  dN  and  [2/SG]  Is  the  fraction  of  grains  in  which  a second 
nucleation  wTll  neighbor  the  first  (with  grains  aligned  in  the  direction  of 
crack  propagation).  We  further  modify  (24)  to  incorporate  the  fact  that  the 
cracks  must  lie  within  a capture  distance  x (Fig.  32)  if  coalescence  is  to 
occur.  Pd( N)  which  results  Is  the  probability  of  nucleation  of  a second  crack 
in  a grain  neighboring  an  already  present  microcrack  in  the  Ntfl  cycle  and  with 
sufficient  proximity  that  coalescence  will  occur  given  sufficient  crack 
propagation. 
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Fig.  32  Maximum  nucleation  site  separation  parallel  to  the  stress 

axis  for  which  coalescence  is  likely  referred  to  as  a capture 
distance,  X. 
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Pd(N)  = P-(N) 


(D) 


X 

IT 


(D)  dD 


(25) 


As  Is  described  in  4,1  Eq.  (25)  and  its  appropriate  modifications  are  used  to 
select  , Dj , W.j  and  which  define  the  nucleatlon  parameters  for  the  second 
nucleatlon.  We  then  take  as  Nmax  for  the  first  nucleation  and  use  the 
single  crack  equations  to  define  the  parameters  for  the  first  nucleation. 


Additional  probability  functions  are  defined  giving  the  probability 
of  nucleation  of  pairs  of  cracks  with  one  Intervening  grain,  etc.  The  form  of 
the  probability  is  the  same  differing  only  In  the  value  of  x,  which  presently 
must  be  determined  experimentally.  Microcrack  growth  of  each  crack  is 
followed  as  per  the  procedure  for  single  cracks,  until  coalescence  occurs  at 
which  point  the  growth  of  the  unified  crack  is  simulated. 

In  Fig.  33  we  show  an  example  of  predicted  to  measured  mean  lifetime 
for  the  two  A1  7075  alloys  discussed  in  Section  4.1.  The  predictions  are  made 
using  the  one  dimensional  simulation  of  nucleation  described  earlier  and  a 
simulation  incorporating  potential  for  coalescence  of  two  microcracks. 

Pd(N)  = 0 for  low  cyclic  strain  and  only  single  nucleation  events  participate 
in  failure,  but  for  larger  cyclic  strain  Pd(N)  / 0 and  coalescence  influences 
fatigue  lifetime.  The  simulated  data  points  in  Fig.  33  are  mean  values  of  100 
simulated  smooth  bar  specimens.  Parameters  used  in  the  nucleatlon  equation 
were  selected  to  fit  the  two  highest  strain  points  for  the  large  grain 
material.  The  predictions  for  the  small  grain  alloy  are  a direct  result  of 
the  simulation  with  no  additional  adjustable  constants. 


Scatter  in  Fatigue  Lifetime 

Scatter  in  fatigue  lifetime  is  determined  in  the  simulation  by 
repeated  running  of  the  simulation  for  the  same  stress  amplitude,  material 
orientation,  etc.  The  scatter  results  from  statistical  variations  in  alloy 
microstructure  from  run  to  run  simulated  by  random  selection  governed  by  the 
grain  and  intermetal  lie  size  distribution  functions.  Our  project  funding 
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concluded  before  an  orderly  investigation  of  the  effect  of  modeling  parameters 
on  lifetime  scatter  could  be  attempted,  so  we  have  no  significant  results  to 
report  in  this  regard.  Presently,  we  jlan  to  defer  a first  test  of  predicted 
to  measured  lifetime  scatter  until  a rare  sophisticated  simulation  based  upon 
the  two-dimensional  nucleation  model,  described  briefly  in  4.1,  is  fully 
operational.  As  an  addendum,  however,  we  want  to  point  out  one  Interesting 
effect  observed  in  the  simulation  of  fatigue  lifetime  scatter  illustrated  by 
the  results  In  Fig.  34.  Lifetime  in  the  case  simulated  is  for  high  cycle 
fatigue  involving  no  crack  coalescence.  A shift  in  the  mean  specimen  lifetime 
is  predicted  with  decreased  specimen  surface  area.  This  "size  effect"  results 
from  a decrease  In  sites  available  for  nucleation,  which  alters  the  statistics 
of  the  failure  process. 


RELATIVE  PROBABILITY  DENSITY 
OF  FAILURE 


5.0  SUMMARY 

A new  approach  to  modeling  of  the  effect  of  alloy  microstructure  and 
environmental  factors  on  mean  and  scatter  in  fatigue  lifetime  has  been 
reported.  The  Monte  Carlo  simulation  of  the  fatigue  failure  processes  de- 
scribed relies  upon  analytical  models  of  microscopic  failure  processes  com- 
prising the  initiation  stage  of  fatigue  failure.  Comparisons  between  pre- 
dicted and  measured  fatigue  behavior  has  been  encouragingly  good  and  we  can 
foresee  that  in  years  to  come  the  general  approach  to  lifetime  modeling 
defined  here  will  have  wide  applications.  The  main  accompl ishments  of  the  re- 
search have  been: 

1)  Test  and  verification  of  a model  relating  alloy  grain  size 
for  7000  and  2000  series  alumni num  alloys  to  the  rate  of 
surface  crack  initiation  at  surface  particles. 

2)  Demonstration  that  crack  closure  stress  is  a principal 
factor  governing  microcrack  growth  in  fatigued  aluminum, 
combined  with  definition  and  verification  of  an  analytical 
model  relating  grain  size  to  closure  stress  and  hence  to 
crack  propagation  rate. 

3)  Definition  of  the  role  of  humidity  in  several  nucleation 
and  microcrack  growth  processes  in  aluminum  alloys. 
Including  development  of  a technique  to  measure  the  effect 
of  humidity  on  the  ductility  of  the  surface. 

4)  Development  and  initial  assessment  of  a Monte  Carlo 
procedure  to  predict  a)  the  mean  and  scatter  in  fatigue 
lifetime  of  certain  aluminum  alloys;  b)  microscopic 
cracking  parameters  of  potential  use  to  evaluation  of 
nondestructive  lifetime  measurement  techniques. 
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Options  to  extend  and  refine  the  models  and  techniques  described  are 
substantial . 

1)  The  nucleatlon  expression  can  be  extended  to  include  the 
effect  of  subsurface  location  or  trl axial  stress  (such  at 
notches)  on  microplastic  ductility  and  hence  cycles  to 
nucleatlon. 

2)  There  Is  the  prospect  that  a systematic  study  of  the  Inter- 
relationships between  alloy  composition,  heat  treatment, 
environmental  humidity  and  surface  ductility  will  lead  to  a 
method  to  predict  the  material  constants  used  In  the 
nucleatlon  equation. 

3)  Slip  distance  models  of  the  type  employed  for  “at  particle" 
nucleatlon  can  be  formulated  to  describe  crystallographic 
and  at  grain  boundary  nucleation,  and  tested. 

4)  With  respect  to  microcrack  growth,  computer  simulation  of 
cyclic  deformation  In  a plastic  zone  confined  to  the  grain 
at  a microcrack  tip  should  provide  substantial  insight  into 
the  Interrelationships  between  alloy  properties,  humidity, 
and  the  material  constants  which  appear  In  the  closure 
stress  equations. 

5)  Methods  to  describe  shear  mode  propagation  such  as  In  the 
presence  of  large  Mode  I closure  stress  must  be  developed. 

6)  Substantial  Investigation  to  develop  models  of  crack  co- 
alescence with  crack  tips  in  common  grains  is  required  and 
should  concentrate  on  the  complex  changes  in  closure  stress 
In  the  region  of  Interaction  which  are  known  to  occur. 
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Comparison  of  Monte  Carlo  generated  coalescence  rates  to 
measured  ones  can  be  used  to  refine  the  description  of  the 
coalescence  process  as  initial  models  are  developed. 

7)  Opportunities  to  test  and  refine  the  Monte  Carlo  simulation 
abound  - considering  that  a more  sophisticated  simulation 
can  Include  two-dimensional  nucleatlon  modeling  Incorpora- 
ting effects  of  grain  orientation  and  alloy  texture,  as 
well  as  orientation  of  the  material  and  stress  axis 
relative  to  the  surface.  By  comparing  predicted  to 
measured  fatigue  properties  on  a microscopic  level  the 
simulation  can  be  used  to  refine  Itself  by  a combination  of 
empirical  means,  and  by  definition  of  areas  requiring 
further  research. 

8)  One  approach  to  the  nondestructive  measurement  of  remaining 
fatigue  lifetime  Is  to  determine  (such  as  by  acoustic 
harmonic  generation)  the  distribution  of  mlcrocraks  on  a 
surface.  The  simulation  can  be  used  to  evaluate  such 
measurement  techniques. 
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